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SUMMARY 


The analysis of a structure composed of beams 
integral with a slab is properly a three dimensional 
blem. In this paper, the need for a three dimensional 
treatment is obviated by an analysis of the slab 
according-to thin plate theory and of the beams, 
treated separately from the slab, by the theory of 
simple bending. The fact that the beams are actually 
integral with the slab is then recognised by enforcing 
displacement and strain compatibility at the slab- 
beam junctions. Both the stress-function, describing 
forces in the plane of the slab, and the transverse 
deflection of the slab are governed by a biharmonic 
equation ; series solutions for these two quantities 
are obtained subject to typical, specified, boundary 
conditions and to conditions of compatibility at the 
slab-beam junctions. The calculation of the coefficients 
in the series solutions would be extremely laborious 
if carried out by hand-computation so that a digital 
computer programme has been written (for the Ferranti 
Pegasus Computer) by use of which the complete 
computation can be carried out in a few minutes ; 
thus the effect of any prescribed variation in the 
parameters describing either the geometry of the 
structure or its loading can be quickly followed through. 

Details of the computer programme are given. 

Further papers will deal with different arrangements 
of slabs and beams, and also with beams which are not 
necessarily of constant depth across the span. Stiffened 
foundation mats will also be considered. 

1, In a compound structure of beams integral with a 
slab, under transverse load, the free bending of the 
beams is to some extent reduced, both by the bending 
of the slab and, when the neutral plane of the slab 
does not contain the neutral axes of the beams, 
by forces induced in the plane of the slab as a result 
of the bending. Some authors!.?2 have suggested that 
this restriction on free-bending of a beam can be 
Simulated by increasing the geometric second- 
moment of area to some ‘ effective ’ value and that, 
for this purpose, it is convenient to think of the beam 
bending, not by itself, but in conjunction with 
flanges taken from the slab. Thus, by concentrating 
attention on the beam-bending, we arrive at the 
concept of ‘effective width’ in which a certain 
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proportion of the slab acts, in bending, with the 
beam. The advantage of this artifice, if we can find 
the magnitude of the effective width, is that design 
curves, which already exist, can be used in designing 
the structure. A difficulty of this approach is that 
not only has the effective width to be calculated, 
but also some basis for its calculation has to be 
established. Such a basis of calculation could be, 
for example, the equality of the maximum deflection 
of the beam, as calculated from some more or less 
rigorous theory, with that calculated from simple 
bending theory using an effective second moment of 
area. An alternative basis could, similarly, be the 
equality of maximum stress. If the effective width 
can be calculated for a few values of the geometric 
parameters describing the structure then design 
curves could be drawn, for a particular loading 
system, from which interpolated widths could be 
read for other values of the geometric parameters. 
It should be noted that the usefulness of the effective 
width concept lies in the assumption that the 
effective width is constant along the span of the 
beam, and of such a magnitude as to satisfy some 
condition, laid down for its calculation, at possibly 
one point only along the span. 

It is probably true to say that the need for such 
a concept has been felt because the amount of 
computation required in the rational design of a 
composite structure was otherwise prohibitive. 
When, however, a digital computer programme is 
available the amount of time taken in computation 
becomes negligible, and an accurate evaluation of 
deflections and stresses can be made without any 
need for the concept ; and, moreover, it becomes 
possible for the designer to modify the size and shape 
of the structure, also the loads which it is to be 
called upon to carry, and to observe, almost 
immediately, what effect his modifications will have 
on deflections and stresses. 

The analysis of a composite slab-beam system is 
really a three-dimensional problem. In this paper 
the three-dimensional nature of the problem is by- 
passed by an initial separation of the slab from the 
beams, analysing the former by thin-plate theory 
and the latter by simple bending theory. The 
structure is then made composite by forcing the 
thin-plate theory and the simple bene rates to 
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slab-beam junctions. Of course, there is some error 
in this process (which will be pointed out in due 
course) but a considerable volume of experiment! 
verifies that the error is not significant. 

In the analysis which is presented here, the 
composite structure, which is supposed to be simply- 
supported at its two edges normal to the spanwise 
direction and to be continuous over several supports 
(Fig. 1), is firstly analysed under the assumption 
that it has no such intermediate supports. With the 
loading specified the computer programme solves the 
problem under this assumption and, in particular, 
produces values of the deflection of the beams at the 
positions where the intermediate supports are in 
reality to be provided. A second run through of the 
programme with unit point loads, only, acting on all 
the beams at one intermediate support (e.g. at 
y = S,L in Fig. 1) again produces values of deflections 
at this and all other support points. This second step 
is repeated with unit point loads acting in turn at as 
many different support positions as exist. A linear 
combination of all the solutions which have thereby 
been obtained can then be made to produce zero 
deflection at all the support points. The complete 
sequence of operations is carried out automatically 
within the computer once the programme and the data 
have been fed into it ; the final solution, printed by the 
machine, is that for the actual structure which is 
continuous over the given number of supports, 
simply-supported on two edges and carrying the 
specified loading. From the method of solution it is 
evident that any postulated sinking of the supports 
could be accounted for, as could a situation where 
the supports are elastic. 

This paper is concerned with the general problem 
of the analysis of multi-span, composite structures 
but there are two particular cases which are of 
special interest. The first is the case when there are 
no intermediate supports; we then have a single- 
span structure, simply-supported at its ends. This 
configuration can be ebtained by giving the s-values 
(describing the positions of intermediate supports, 
Fig. 1) appropriate magnitudes. The second interest- 
ing particular case is obtained by considering a 
structure with a number of spans, say five, which is 
loadcd symmetrically with respect to the inter- 
mcdiaic supporis; the centre span of the five then 
effectively behaves as a single-span structure with 
its beams encastré at the edges. The slope of the 
slab, in the spanwise direction, is also zero at the 
edges, although its deflection is not forced to be zero 
there by the above considerations, so that we do not 
by this artifice quite get the solution to the problem 
where the whole structure (slab and beams) has 
encastré edges; the beams do, however, gencrally 
supply the greatest part of the resistance to bending 
inthe spanwise direction and the fact that they can 
be made encastré means that a good approximation 
is achieved to the problem where the whole structure 
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has encastré edges, the approximation becoming 
better as the beam spacing is reduced. To analyse a 
structure falling within either of the above special 
cases no modification to the computer programme 
needs to be made ; it is merely necessary to arrange 
for the data-tape to contain appropriate essential 
information. 


2. Fig. 1 shows the arrangement of the structure and 


the meaning of the symbols used. It is assumed that 
simple-support conditions exist on y = 0 and y = L, 
The beams in the structure are of rectangular cross- 
section, are equally spaced and, as far as the analysis 
és concerned, are infinite in number. That is to say, 
the analysis is strictly valid only for those parts of the 
structure which are uninfluenced by its termination 
in the x-direction. Three kinds of loading have been 
considered in writing the programme, 

(i) Uniformly distributed load on the slab (dead- 
weight or live-load), 

(ii) Uniformly distributed load along the beams, of 
the same intensity on each beam (dead-weight 
or live-load), 

(iii) Groups of point loads on all the beams, of the 
same value on each beam at any one spanwise 
position. (It is necessary to consider several such 
groups of point-loading). 

With these kinds of loading, model tests! have 
shown that beams two spaces from the termination 
in the x-direction behave as though the beams were 
repeated to infinity, and comparison with the model 
test results shows that our analysis is also within 
practical accuracy for beams one space from the 
termination. 


3. The equations governing the plane stress-function, 


f, and the transverse deflection, w, are respectively 
Vif=0 
and ° ° ° ° ° ° (1) 
where D is the flexural rigidity of the slab and ¢ 
the intensity of uniformly distributed load acting 
on it. Solutions to the equations (1) in the region 
—al2<¢x <¢ a/2,0 < y < L, between two adjacent 
beams, may be taken in the form 
f = <Xsiney[Ancoshex+ Byxsinhex+Cypsinhex ) 
+ Dyxcoshex), | (2) 

w= Lsiney[gn+ Pacoshex+ Opnrsinhex 4 

+ Rysinhex+Spxcoshex), J 
the summations extending over positive integral 
values of m. By virtue of the symmetry about both 
x = 0 and y = L/2, these reduce to 
f= Lsiney[ A pcoshex+ Boxsinhex], \ Z (3) 
w= Lsiney[(gn+ Pacosher+Onxsinhex), 
where ¢ = mn/L, gn = 4qL4/(nx)5D, and the sum- 
mations extend over odd positive integers m. The 
series (3) satisfy the necessary conditions of symmetry 
and also three of the four conditions of simple 
support on y = Oand y = L, viz., the conditions that 
on these edges 

(i) the deflection, w = 0, 
aw aw 


(ii) the bending moment, My = Px +o Be = 0, 


(iii) the normal stress in the plane of the slab, 


The fourth condition, which is not however satisfied, 

is that 

(iv) the shear stress in the plane of the slab, 
an arf 

x by 


= 0, 












wa 
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At the expense of considerable complication in the 
first of (3), the series form for the stress-function, 
this fourth edge condition can be satisfied and the 
solution made mathematically correct, but model 
tests! have shown that theoretical results based on 
(3) are, within the limits of experimental accuracy, 
in reasonable agreement with experimental deflection 
and strain measurements. Theoretical work using 
a series solution (other than the first of (3)) which 
does allow satisfaction of condition (iv) has also 
shown that little error results from using the first 
of (3) for the stress-function*. 

*The question of the satisfaction of condition (iv) is discussed in 

an Appendix to this paper, where it is demonstrated that the 

failure of the series (3) formally to satisfy the condition is not 
significant, 

4. The constants An, Bp, Pp and Qp in the series (3) 
remain to be found by satisfaction of further condi- 
tions at the slab-beam junction. First, conditions of 
symmetry there require that 


Le 
ox 
= on x = 4/2, forall y, - - (4) 
u=0 


where # is the displacement component in the x- 

direction. Satisfaction of the first of (4) by (3) leads to 
1 a ea 

To satisfy the second of (4) we employ the relation 


E = =xx—oyy, and the fact that « is zero on 


x = 0, for all y, to obtain 


I—o 1 @ ea 
A, = Ba [ize 3-3 | (6) 














Two further relations are needed before all the 
constants can be computed, and these are obtained 
by consideration of conditions of equilibrium and 
of strain-compatibility at the slab-beam junction. 
By reference to Fig. 2, consideration of the equili- 
brium of a beam element (dy) gives 
pidy + dF —2S,;dy =0, 
aN — 2hxydy =0 } (7) 
and @M+ Fdy— 2hxyldy = 0, 
where J] = (4 —ny/2 M = El a 
’ dy? 
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The quantity #1, tne transverse load per unit iength 
acting on the beam element, includes both externally 
applied loads and the deadweight of the beam. 
From the first and third of (7), by elimination of 
F, we obiain 
day 
oy 
Equality of strain at the slab-beam junction gives 
(it being assumed that the direct force N is distributed 
uniformly over the beam cross-section) 

N aw — — 

ta + El oa =I — oe. ° ° - (9) 
The second of (7), together with the fact that N = 0 
at y = 0, enables us on integration to replace N in 


(9) by —2n . Satisfaction of (8) and (9) by (3) 


otw 
EI dyt + 2S; --p! — 2hl => 0. ° (8) 


then requires that 
P,[ EI e*?+2D(1—e) T]+ Qn[4$EI 2a T—2D(1+0)T 
+ Dea(1—o)]—Ap([2hleT]—B,[2A1T + hlea] 
Ka 4gEI 4p 


Cet LDCe8 * EC’ as 


and 
Pal El] + Qnl}ElaT) + 4al a + (14 °)| 


MT , ah , 2 ; 4g El 
+ Bal Seg t Sat i+ il+oer| ~~ Cpe 


. . (11) 
where T = tanh}ea, C = coshhea, p is the uniformly 
distributed part of #1 and Ky is the Fourier cofficient 
in the expansion of the remaining part of #1 in odd 
sines. In order to find Pp, Qn, An and By for any 
given system of loading (symmetrical about x = 0 
and y = L/2) it is necessary to solve (5), (6), (10) and 
(11) simultaneously for as many successive values 
of m as may be necessary to obtain practical con- 
vergence of the series (3). This kind of repetitive 
calculation is ideally suited to a digital computer 
and a previous programme has been written, and 
described elsewhere*, for composite slab-beam 
structures spanning between simple supports. 





‘ 
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Fig. 3. 


5. When the composite slab-beam structure has 


intermediate supports across the span, the previous 
analysis and computer programme would still 
suffice if the values of the redundant reactions at the 
support were known ; for these reactions could then 
be taken into account in the terms Ky. These 
reactions are not known, however, but instead must 
be found as part of the solution. As an illustration 
of the computational procedure, incorporated in a 
new programme, consider the arrangement, shown 
in Fig. 3, in which there are four intermediate 
supports. In addition to uniformly distributed 
load over the slab (of area-intensity g) and uniformly 
distributed load along each beam (of line-intensity p) 











we make provision for point-loads, Wo, W; and Wg, 
acting on each beam at the centre-points of the 
individual spans. The redundant reactions, also 
shown in Fig. 3, are treated as negative loads and 
denoted by — Ws and —W,. If we consider a 
single point-load, W, acting on each beam of a 
simply-supported structure of span L at a distance 
sL (0 < s < 1), from the edge y = 0 (Fig. 4) it may 
be shown that the coefficients Ky are equal to 
(2W/L) sinnsx. When all the point-loads and 
redundant reactions of Fig. 3 are acting the first 
term on the right-hand side of (10) becomes 
2[W osingnn+ Wy | singn(si+se)x-+-singn(2—s;—sq) | 

+ Wea sindnsyx+sinjn(2—s})x } 

+ Ws| sin nsyx+sin n(1—s,)7} 

+Wa4| sin nsax+sin n(1—sg)n}]/LCe® + (12) 
in which, of course, the arithmetic values of Ws and 
Ws, if we knew them, would be negative. If we 
remove the intermediate supports, that is put 
Ws = W4= 0, a first run-through of the computer 
programme would provide values of the transverse 
deflections at y = s,Z and y = seZ under the action 
of the loads g, ~, Wo, Wi and We; let these 
deflections be d and D, respectively. A second run- 
through of the programme with g = ~ = Wo = W, 
= We = W, = 0but Wz = 1 will give corresponding 
deflections (d’ and D’) due to unit, negative, support 
reactions at y=s,l. A third run-through with 
gq=p = Wo= W, = We = Wz = 0 but Wa = 1 
will give corresponding deflections (d’”’ and D’’) due to 
a unit, negative, support reaction at y = sel. If we 
now define « and 6 by 


d = ad’ + pa” 


and run the computer programme through a fourth 
time with g, ~, Wo, Wi and Wg taking their given 
values, W3 = — a and Wg = — 8, the computer 
will obtain the solution with zero deflection at both 
y = s,Landy = sel and, on account of the symmetry 
about y = L/2 which has bee1u preserved throughout, 
at the other two supports as well. It is clear that 
any specified sinking of the supports could be 
accounted for in the foregoing method, as could a 
problem with elastic supports. The computer pro- 
gramme is so designed that the four runs-through 
described above follow automatically, in succession, 
after a single, original, input of t. e . ecessary data. 
The intermediate results, d, D,d’' D d” and D” are 
retained in the computer which also automatically 
solves equations (13) to obtain values for « and 8. 


6. Working details of the computer programme will 
be published elsewhere; it is sufficient here to 
make a few general remarks on the input of informa- 
tion about the structure and its loading, and on the 
output of results. 

The input tape, for four intermediate supports, 
consists of the sixteen, signed, numbers 

a/2, b, d, L, h, E, SG, q; Wo, Wi, Wa, p, $1, S2, Ws, W4. 
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The last two numbers are not really part of the input 
information and are always put equal to zero, their 
values for the first run-through ; they are included 
in the data-tape for subsequent ease in the organi- 
sation of the computation. 

The output tape is printed as follows 


Q1 


Pi 
wat y=L/2, for n=1 
Qs 
for ~,¢7, Wo, Wi and WV, 


having their given values, 
and Ws = Ws ed 0, 


Ps 
wat y=L/2, for n=3 


Qis 
Pis 
waty=L/2form=15 J 
the twenty-four above quantities for 

p=q=W o=W1=W2=W4=0, Ws=1, 
the twenty-four above quantities for 

p=q=W o=W1=W2=W3=0, Wa=1, 
the twenty-four above quantities for 4, Q, Wo, Wi, 
W3 as given, and Ws=—«, W4= —8, 





waty = 1 
wat y = sel 
The last two quantities printed should, of course, 
be zero and they are included in the output as a 
check that the computation has been carried out 
correctly. The computation could be carried out 
up to any value of but exploratory work has shown 
that, for the types of loading considered, eight 
terms in the series solutions is quite enough. The 
coefficients Q, and Py, are printed so that, with a 
small amount of hand computation the transverse 
deflection can be calculated from (3) at any point in 
the structure. The transverse deflection is calculated 
by the programme for y = L/2 and printed so that 
if it is desired to analyse a single-span composite 
structure of this kind, for which the beam-ends are 
encastré, the present programme (without any 
modification) will produce the centre-span deflection. 
All that need be done to achieve this effect is to put 
Wy = We = Wo, 1 = 0-2 and sg = 0-4 on the 
data-tape, with # and g equal to their specified 
values. We shall then obtain the centre-span 
deflection for an encastré structure of span L/5 under 
uniformly distributed load on the slab, uniformly 
distributed load along the beams and a point load 
(W 9) at the centre-span of all beams. 


for the last loading condition. 


. As an illustration of the use of the programme we 


consider a single-span composite structure and 
investigate the difference, in the centre-span 
deflection, between conditions of simple support and 
conditions of fixity of the beam-ends. For this 
purpose we consider two different specimens whose 
geometry is defined as follows : 





a/2 





Specimen 1 3-625 








Specimen 2 | 6-0 





























All dimensions are in inches and E in tons/in®. For 
Specimen 1 loading conditions are taken as p = g = 0 
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with W, the point-load at the centre-span position 
of each beam, as 60 lb. For Specimen 2, W = 0, 
p = 0-723 Ib/in. and g = 1-922 Ib/in?. Such speci- 
mens, with the loads quoted, have been studied 
experimentally! and measurements made of the 
mid-span deflections (wa) of the centre beam of 
such structures with five beams composite with the 
slab when the edges were simply supported. In these 
experiments the model structure was built up from 
aluminium sheet and bar, and the beams cemented 
to the slab with araldite resin. From subsidiary 
experiments it was found that this method of 
construction did reproduce the conditions assumed 
in the theoretical work up to the intensity of loading 
to be used. The point-loads, which were applied 
to the centre-span point of each beam, were trans- 
mitted by a system of levers from dead-weight 
loading pans located on either side of the model. 
The uniformly distributed load was produced by 
compressing an air-bag between the model and a 
thick steel plate. The air-bag could not, in fact, 
be made to cover the very corners of the model 
without overhanging in other parts but that this 
system did produce an effective uniformly distributed 
load was proved by applying it first to a thin slab 
with stiffening beams and comparing experimental 
deflections with those which had been calculated 
from accepted theories. Calculations have also been 
made, using the earlier computer programme? to 
obtain those deflections (wp), again for the case of 
simply supported edges. We now use the computer 
programme which is the subject of this paper to 
calculate the mid-span deflection (we), when condi- 
tions of end-fixity apply, for the same specimens and 
applied loads. In order to do this, the data-tape 
must be punched with the information in the table* 
following, the first line being for Specimen 1 and 
the second line for Specimen 2. 


*In this table the loads are quoted in ton units. 


L h E c 
145-0 | 0-15/4777-0 | 0-3 
0-9 120-0 | 0-15/4777-0 | 0-3 
We p s1 | Sg Ws 

0 0-2 0-4 0 


0-2) 0-4 0 


The results are collected and compared in the 
following table, the entries being in inches x 10-%: 





Wa Wb We 





Specimen 1 51 47 14 





Specimen 2 66 64 38 




















The small discrepancy between wa and wp for 
specimen 1 can be attributed, in some degree, to 
experimental error and also in part to assumptions 
made in the theoretical work (such as the assumption 
that the direct force N is uniformly distributed 
across the beam cross-section). The effect of fixing 
the beam-ends against rotation is to produce a 
considerable reduction in the value of the deflection 
although, for these specimens, the beam spacing 1s 
large in relation to the other dimensions. 
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8. The computer programme described in this paper 

deals with four intermediate supports, symmetrically 
disposed about y = L/2, but it can be immediately 
applied, by suitable punching on the data-tape, to 
problems where there are no intermediate supports 
at all. If there is an even number of intermediate 
supports greater than four, or if there is an odd 
number of intermediate supports, again symmetrically 
disposed about y = L/2, a small modification to 
the programme as at present written would be 
necessary. 
The time taken for computation on one set of data 
is approximately three minutes, which means that 
possible alternative values of parameters can be 
followed through in a very short time. 


~ 
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APPENDIX 


The series form for the stress-function, the first of (3), 
has been shown to satisfy all the necessary conditions 
except the one requiring %j to be zero on the simply- 
supported edges of the slab. The analysis in the body 
of the paper is mathematically lacking to this extent, 
but it has been assumed that no significant error would 
arise on this account ; the validity of this assumption 
is supported by the agreement between our results and 
those of a considerable body of experiment!. We now 
derive an alternative, more complicated, series-form 
for the stress-function, f, in order that a theoretical 
test of the assumption can be imposed. It is convenient 
to change the origin of coordinates (Fig. 1) to the 
point x = 0, y = L/2. It can be then verified that the 
series 

f = Xy[Arcoshex+ Byxsinhex]}cosey 

+ XpCp[4lsinhtaLcoshay—ysinhaycoshtaL]cosax 

. . . (14) 
where ¢ = rx/L, « = nz/a, and the summations extend 
over odd integral values of r and m respectively, 
satisfies all the conditions which (3) had satisfied and, 
by virtue of the introduction of the Cy terms, also 
allows xy to be made zero on the simply-supported 
edges, y = +L/2. The corresponding series form for 
the transverse deflection is 

w= Xlgr + Preoshex + QOrxsinhex]cosey 
where gr = 4gL4sindrx/(rx)5D. 

For xy to be zero on y = +L/2, for all x, we obtain 
YnCnaltal + coshtaLsinhhaL)sinax = 

Xr! Bre(sinhex + excoshex) + Arye*sinhex)|sinkrnx - (15) 
In order to find the coefficients Cy from (15) it is 
necessary to expand both sinhex and xcoshex, on the 
right hand side, in odd sine series (of sin ax). If the 
double summation on the right-hand side of (15) is 
then re-arranged we obtain equations of the form 

Cr = UrbmnAr + UrymPBr 
for each value of n. 

When the C, terms were absent the four coefficients 
Ay, By, Py and Q;y could be calculated for ry = 1 and 
the calculation repeated in turn for r = 3, 5,7,..., for 
as many values of 7 as were found necessary, each 
group of four constants being calculated separately ; 
in fact, we had a system of four simultaneous equations 
to solve for each value of r. The effect of the inclusion 
























of the C, terms is to make all the constants in:er- 
dependent. Thus, if we wish to calculate k terms of 
the series, we have to solve, not & independent sets 
each of four equations, but instead one set of 5k 
simultaneous equations. The itude of such a 
computation using hand calculating machines is 
formidable but for reasonable values of & it presents 
no difficulty when a digital computer is employed. 

In the absence of the Cy terms the equations (5), 
(6), (10) and (11) expressed satisfaction of various 
conditions for all values of y and it was from these 
four equations that we were able to calculate the 
values of the constants. The inclusion of the Cy 
terms modifies all but one of these equations by the 
inclusion of functions of yy such as coshay and ysinhay. 
Before the equations can now be satisfied for all y it is 
necessary to expand coshay and ysinhay in odd cosine 
series (of cosey). For example, the corresponding 
form of (6) is 

Er[{ Br(1—o)—A re(1+-0) }sinh}ea 

—4(1+0) Breacosh}ea]cosey 

+ EnCnsingnn|| $(1-+0)aLsinhgaL 

—2coshtaL \coshay—| (1+) acoshtal \ysinhay] = 0, 


(16) 
and we find that 
coshay = = Ph. ee sin = cosh =. cos 
") x(n8L2 aes) "2D 2 "y 
oe ee 
2Lra2 +. or 
y sinhay = =r Ena sine sinh > 
4Lan h aL 
= (m2? a%2) cos. z cosey. 
(18) 
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If we let the quantities in curled brackets, in (17) and 
(18) respectively, be represented by J and JJ then (16) 
becomes 


Ap{—esinh}ea] + By liz sinh tea — tea cosh bea] 
2 


+ InCe sin dre (de sink al — >— cosh tal 1 


— { acosh tal | n] = 0. 





The equations corresponding to (10) and (11) can be 
obtained in a similar manner. Equation (5) is un- 
affected by the inclusion of the Cy, terms. 

The experimental and theoretical work mentioned 
in the body of the paper suggests that the inclusion of 
the Cy, terms will make no appreciable difference to 
the results obtained and that, because of this, the simple 
theory can be used with confidence. To check this 
suggestion a digital computer programme has been 
written which includes the Cy terms and the computa- 
tion has been carried out on data for which the solution 
by the simpler theory had previously been obtained by 
another programme’. The problem is that of a single- 
span composite structure of a slab and beams simply 
supported along the two edges normal to the spanwise 
direction and loaded equally at the centre-span points 
of all beams. The coefficients which were obtained with, 
and without, C, terms are collected in the table below: 
It is clear that the inclusion of the Cy, terms does not 
make any practically significant difference to the 
values of the coefficients. 


Discussion 

The Council would be glad to consider the publication of 
correspondence in connection with the above paper. Communi- 
cations on this subject intended for publication should be 
forwarded to reach the Institution by August, 31st. 1961. 
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Wind Structure and the Assessment 
of Wind Loads 


by Walter A. Woeber, D.Sc., M.I.Mech.E., M.(S.A.)I.C.E. 


1. Introduction 


Structural analysis devolves upon two essential 
operations : the assessment of the external forces to 
which a structure may be subjected, and the deter- 
mination of the internal stresses induced by these 
forces in the various portions of a structure. 

In specifying the external forces, recourse is usually 
had to accepted codes of practice which recommend 
the loading to be taken into account for a variety of 
conditions. It is not sufficiently realised, however, 
that these load classifications are of necessity imperfect 
representations of actual conditions. This arbi- 
trariness is particularly evident in the case of the 
design wind loads laid down in building codes. 

The purpose of this paper is to bring the problems 
associated with the assessment of wind loads in the 
design of exposed engineering structures within the 
purview of that part of dynamical meteorology which 
deals with atmospheric motion near the ground. 


2. Atmospheric Turbulence 


Atmospheric motion, which becomes apparent to 
the observer as wind, is generally of a very complicated 
nature. The term ‘wind’ is restricted to the 
horizontal component of this motion and is specified 
by its speed and the direction from which it blows. 
Wind speed and direction are subject to rather wide 
and rapid fluctuations and exact specification or 
mathematical analysis is difficult. 

The motion of the atmosphere, particularly in the 
surface layers, is usually turbulent. Firstly, there is 
the mechanical effect of irregularities on the ground, 
giving rise to frictional eddies which are carried along 
with the wind. Secondly, owing to the heating of the 
air at the earth’s surface, convectional eddies are 
set up and develop freely when the heating is intense. 

Dry air, when heated by contact with a warmer 
ground surface, rises and expands and cools adia- 
batically at the rate of 1°C in 100m. This is known 
as the adiabatic lapse rate. 

By day, the atmospheric temperature gradient 
may exceed the adiabatic lapse rate and facilitate the 
formation of turbulence. The thermal currents due to 
convection may either form part of the general 
turbulence or appear as small travelling eddies caused 
by local temperature differences on the ground, where 
bare soil, rock or sand absorb the incoming solar 
Tadiation more readily than grasslands, wooded areas 
or water surfaces. 

By night, when the ground is cooled by the flow 
of the outgoing radiation, the air tends to attain the 
Same temperature as the ground with which it is in 
contact and the lapse rate may become inverted or 
isothermal through a layer of 500 ft. or more. An 
inversion is characterised by its pronounced stability 
Which damps out convection and prevents turbulent 
mixing. 


The degree of turbulence which may be present 
in the atmosphere at any time is clearly demonstrated 
by the trail of smoke from a chimney. Disregarding 
for the purpose of this discussion, the force of 
buoyancy of hot plumes and the topographic effect 
on plume dispersion, the shapes traced out in a fresh 
breeze will assume distinct patterns which can readily 
be associated with the various temperature lapse 
rates, as shown in Fig. 1. 
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Fig. 1. 


The variation of atmospheric temperature with 
height is shown diagrammatically on the left hand side 
of Fig. 1 for each of the five cases (a) to (e) illustrated. 
The adiabatic lapse rate has also been plotted to serve 
as a comparator. 
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In 1 (a) the lapse rate is clearly super-adiabatic 
and the marked thermal turbulence resulting there- 
from is demonstrated by the ‘looping’ of the smoke 
emission. 

Temperature conditions approximating to the dry- 

adiabatic lapse rates are represented in 1 (b), and the 
resulting smoke plume is in the form of a horizontal 
cone. 
A temperature inversion in the atmosphere, 1 (c), 
denotes thermal stability, illustrated by the almost 
laminar flow from the stack which is referred to as 
* fanning.’ 

Conditions associated with an inversion near the 
ground and a super-adiabatic gradient above it, are 
shown in 1 (d). The thinning out owing to convection 
from the top of the otherwise laminar emission has 
been termed ‘ lofting.’ 

Finally, in 1 (e), super-adiabatic conditions exist 
near the ground and an inversion above the top of 
the stack. Rising of the smoke is prevented by the 
stable conditions at the upper level; the smoke 
stratifies and forms a screen Ci which atmospheric 
turbulence causes ‘fuming’ or the spreading down- 
wards of the smoke. This is the condition responsible 
for the formation of what has popularly been 
termed ‘ smog.’ 

Steady or streamline motion is a rather exceptional 
case in nature and even the almost laminar emission 
shown in Fig. 1 (c)'appears as a steady motion only 
to the human observer who looks at the plume from a 
far distance. In the case of the laminar emission, the 
radius of the circular cross-section containing 90 
per cent of the suspended smoke particles will be about 
¢in. at a distance of 30ft. from the point source 
(top of stack), and about Sin. at a distance 3,000 ft. 
further downwind. In an unsteady wind the turbulent 
exchange is of a much higher order and the smoke 
spreads out very much faster, as shown in the 
illustrations. 

The transition from a steady to an unsteady wind 
motion has its parallel in fluid flow through pipes, 
well known to engineers, where a change from laminar 
to turbulent flow depends upon the ratio of the product 
of the radius of the pipe, 7, and the average velocity 
ef the fluid, wu, to a physical property of the fluid, 
the kinematic viscosity, v. The ratio, «/v, is known 
as the Reynolds number (Re). 

Remembering that the kinematic viscosity is 
defined as the ratio of the absolute viscosity, u, to the 
mass density of the fluid, p, the importance of the 
Reynolds number becomes immediately apparent from 
the following relationship : 


, dynamic force 
2 a ° ; : 
url/v o $pu*/(uu/r) o viscous force ? 


At low values of the Reynolds number the flow is 
controlled by viscous forces, disturbances are damped 
out and the motion is steady. When this ratio 
increases, the dynamic forces become more powerful, 
disturbances are not so easily damped out and the 
flow breaks down into a violent eddying. 

The analogy of pipe flow is, however, relatively 
simple when compared with atmospheric motion and 
it is not possible to apply the Reynolds number in a 
straightforward manner to problems connected with 
the flow of the natural gusty wind. 

When considering atmospheric motion, the dimension 
r is presumably to be taken to measure the height 
of the homugeneous air mass at the observing station, 
so that motion on such a vast scale must be turbulent 
even at relatively low wind speeds. Only during an 









The Structural Engineer 





inversion or with air motion corresponding to light 
winds up to force 3 on the Beaufort scale (5 m/sec) 
is the flow nearly laminar. 

The above considerations show clearly that the state 
of the atmospheric motion is not specified by the 
Reynolds number alone but that account must also 
be taken of the prevailing temperature gradient, 
Conditions of zero lapse rate, as they occur in a wind 
tunnel, make the interpretation of laboratory tests a 
matter of considerable difficulty. The wind speed and 
viscosity are the same in the wind tunnel and in the 
atmosphere, so that the Reynolds number in the test 
is reduced by the scale of the model. In practice there 
is, however, a vast difference between the two Reynolds 
numbers which clearly indicates that the simulated 
flow in the wind tunnel is far from being dynamically 
similar to the atmospheric motion. 

A criterion for deciding whether steady or unsteady 
conditions exist in the atmosphere and, further, 
whether atmospheric turbulence will subside or 
increase, has been enunciated by L. F. Richardson 
whose work in this respect is comparable to that of 
Osborne Reynolds. For the purpose of this paper the 
Richardson number (Ri) may be stated as follows: 


ni-£(%+1)/(#) , 


where « = average wind velocity, 
z = height above the ground, 
g = acceleration of gravity, 
+ = temperature at level z, 
_ = adiabatic lapse rate, 
dz/dz = observed temperature gradient, 
du/dz = velocity gradient. 

The critical value of this parameter, below which 
turbulence sets in, and above which conditions tend to 
become vertically stable, has been estimated by 
Richardson at unity. Although other investigations 
have thrown some doubt on the validity of Réerit = 1, 
work carried out by C. S. Durst? in the lower layers of 
the atmosphere has shown that there is ample evidence 
to support Richardson’s criterion. 


3. Wind Structure and Gustiness 


From the foregoing it is clear that the motion of the 
atmosphere is not constant, neither in space nor in 
time; wind speed and direction vary continually. 
If these variations take place at definite intervals and 
reach definite amplitudes, they are referred to as gusts. 

The gustiness of the wind is defined by the passage 
of a quantity of air (M) across a control surface (A) : 


wa) pucosads - . . - (3) 
A 


where p = air density, 
“ = wind velocity, 
« = angle of incidence of wind direction, 


f do = A = area of the control surface. 


Depending upon whether u, a, p, or A varies, the 
above expression permits distinctions to be made 
between sudden fluctuations of wind speed, direction, 
density or the cross-sectional area of the wind stream. 
In the present paper attention will be focused on the 
speed and directional variations of the wind as 
recorded by a pressure-tube anemometer. 
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The pressure-tube anemometer was designed by 
W. H. Dines, F.R.S., when, after the Tay Bridge 
disaster in 1879, special attention was devoted to the 
measurement of gusts. The instrument is based on 
the principle of the Pitot-tube and is kept pointed 
into the wind by means of a vane. The static tube is 
vertical. The wind blowing into the impact tube 
produces an increase of pressure inside it which is 
transmitted through the central supporting tube to 
the recording apparatus. The central tube is surrounded 
by an outer tube in which there are four rows of holes 
to provide a suction effect. The difference between the 
pressure and the suction effect is made to actuate a 
float or diaphragm whose motion is recorded by a pen 
on a chart. 

The trace of the pressure-tube anemometer affords 
information as to the internal structure of the wind. 
The width of the trace is an indication of the degree 
of gustiness, and it is usually possible to insert, by eye, 
a mean line through such a trace and so to determine 
the mean wind at any time. 

A closer examination of an anemometer trace reveals 
that periods of stronger and weaker fluctuations 
alternate so that the record assumes the character of a 
stepped profile. This is more clearly brought out in 
the enlarged reproduction of a recorded wind gust, 
shown in Fig. 2, which was obtained by increasing the 
speed of the drum recorder. 
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Fig. 2. 


The short fluctuations, with periods of a few seconds 
only, are caused by surface friction ; superimposed on 
these and characterised by comparatively wide fluctua- 
tions are disturbances with periods of several minutes. 
The latter are caused by thermal eddies and are 
largely responsible for the gusts and lulls in the wind 
motion. 

The basic character of the record of a gust is nearly 
always the same: after a steep rise, the wind speed 
dies down in gradual steps, the gust is followed by a 
lul!. Although it might appear that a mean line drawn 
through the succession of irregularities follows a sine 
oscillation with a decreasing amplitude, the numerous 
errors in geometrical form preclude the evaluation of 
anemograms by the methods of analytical geometry. 
It is however possible to apply the impulse theorem 
to the irregular atmospheric motion. 

As stated above, the mean line through an anemo- 
Meter trace can usually be determined with sufficient 
accuracy by eye. For a relatively short time interval, 
Tepresented by the ‘sampling length’ of the mean 





157 





line, the mean velocity can be stated for any point on 
the curve. If the mean values of the components of 
the wind motion, referred to the rectangular co- 
ordinate system in space, are denoted by wu, v, and wv, 
and the temporary deviations from the mean by 
u', v', and w’, then the actual wind velocities at any 
particular time are 
u=—ut+tu,v=v+0,w=w+w. 

In applying the impulse principle to wind motion, 
it is not sufficient to apply it to the average motion ; 
the mean value must be obtained from the individual 
impulses passing through the fixed control surface (A), 
referred to earlier. In what follows, the u-component 
only of the wind motion will be considered. 

During the time interval dt, a mass p udt passes 
through the control surface with the velecity wu. The 
impulse for the time dt is therefore p u2dt. The total 
effect over a longer period of Time T is 


T 
Lf e “2 dt per second. 
) 


Using the notation previously adopted, namely placing 
a bar above the symbols to denote mean values, the 
above expression may be written in the form op u?. 
Now w2 = (wu + wu’)? = (uJ? +2uwu'+(u')2 + (4) 
In forming the mean value of u?, written u?, it must 
be borne in mind that w already represents the mean 
value of «; further, the mean value of w’ is zero by 
definition, but this does not apply to mean values of 
squares of u’ Hence 

put=p(u)i+ oie)? - - -<- = 9) 
Thus, the total effect is obtained by adding to the 
impulse of the average motion the impulse resulting 
from the mean value of the square of the deviations. 


4. Evaluation of Anemograms 


The force on an object submerged in a fluid stream 
depends upon the shape, the size, the velocity of the 
flow, the density of the fluid, and the viscosity. 

By the application of the principles of dimensional 
analysis, the relationship may be expressed as follows : 

F=}tewRo - ; , ‘ - (6) 
where F = fluid force, 
9 = mass density of the fluid, 
“ = velocity of flow, 
1 = typical dimension of the submerged object, 
@ = a function of the Reynolds number. 

The factor 4 has been inserted in equation (6) 
because $9? has a physical significance as the pressure 
which would result from a complete conversion of the 
kinetic energy of the wind into its pressure equivalent. 
Since, in the case of wind motion, the Reynolds number 
is largely influenced by the thermal stratification 
of the atmosphere, 9 is, in a sense, a function of a 
function. The degree of gustiness present in the 
atmosphere at any time is, however, reflected by the 
width of the anemometer trace so that, in the evaluation 
of wind records, due cognisance is taken of Reynolds’ 
and Richardson’s parameters. 

As / is a linear dimension which defines the shape 
of the submerged object, /2 may be replaced by 4, 
its projected area on a plane perpendicular to the fluid 
force. In order to allow for the different resistances to 
the air flow which arise from the particular shapes 
of the objects exposed to the wind stream, the term A, 
for the projected area, is grouped with a shape factor, 
S. This coefficient has widely different values for 
bodies of different shape and is lowest for streamlined 


objects. 
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BLOEMFONTEIN 
ANEMOMETER 
31-8-41 


In wind pressure calculations, where the shape of a 
structure exposed to the wind is unusual and available 
data are insufficient, recourse must be had to wind 
tunnel experiments to assign a value to the shape 
coefficient. The variation of the shape factor with 
the Reynolds number, the so-called ‘scale effect,’ is 
fortunately not appreciable in the case of engineering 
structures, which generally have sharp edges or corners 
so that the shape coefficients obtained under simulated 
flow conditions may be accepted, albeit with caution 
as previously emphasised. 

By taking p to be 0-00238, its value at 15°C and 
760 mm. Hg., and expressing the wind velocity in 
miles per hour, U, equation (6) may be written in the 
following form : 


F=0-000SAU2 - - + + 


In order to allow for the racking action of the 
pulsating wind force, the design wind velocity, U, 
is grouped with a gustiness coefficient, G, so that the 
wind pressure formula to be applied in structural 
design may be stated as follows : 


P=0-0020SGU2 - - - = () 


where P = F/A = wind pressure intensity in Ib/ft.?. 

The gustiness coefficient plays a somewhat analogous 
part to that which the coefficient of impact plays in 
structural design to denote the number which the 
computed static stress is multiplied to obtain the 
value of the increased stress assumed to be caused 
by the suddenness of application of the live load. 

The degree of gustiness is indicated by the width 
of the trace of a recording pressure-tube anemometer, 
as shown in Fig. 3, which has been chosen to illustrate 
two cases of extreme gustiness recorded at the 
Bloemfontein and Port Elizabeth aerological stations. 

S.P. Wing? recommends an increase of 25 per cent 
of the readings taken with a cup anemometer. 





Fig. 3. 
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PORT ELIZABETH 
ANEMOMETER 
26-4-47 


W. A. Woeber‘ relates the gustiness coefficient to the 
ratio of the mean amplitude of the valleys and peaks 
traced by the recording pen to the highest mean wind 
over a specified sampling length of an anemogram. 
Thus 

G=1+4'/u ° ° ° ; - (9) 


It now remains to determine U, the wind 
velocity to be used in structural design for any 
particular locality. 

An examination of anemometer records at any 

observing station reveals that high winds occur 
infrequently ; it is essential, therefore, to use as long 
a period as possible in selecting an anemogram of 
extreme wind conditions for the purpose of evaluation. 
The record so selected should in fact represent the 
wind extremes over a period of not less than twenty 
years. ; 
Exhaustive investigations by W. Schmid® have 
shown that, strictly speaking, wind records are valid 
only for the immediate neighbourhood of the aerological 
station at which they were taken, and that there is no 
ready method of obtaining representative data for 
larger districts from the extreme values so determined. 
It is therefore necessary to analyse the wind records 
and to specify the maximum intensities for a number 
of well distributed aerological stations and thereby 
divide a district into regions according to theif 
respective windiness. A survey of this nature has been 
carried out by W. A. Woeber® for the Union of South 
Africa. 

Since the anemograms selected for analysis record 
exclusively the horizontal component of the wind 
motion, the following considerations will be confined 
to the #-component of this motion. It is clear that the 
total wind effect which must be accounted for in 
structural design is equal to the sum of the effects 
due to the mean wind and the fluctuations about the 
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DURBAN ANEMOMETER 
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mean. In accepting this argument, equation (5) may 
be re-written in the following form : 
put == p US = olut-+w) - - + (10) 

Thus, the two principal quantities to be obtained 
from the anemogram are wu, the ‘highest mean,’ 
and #’, the mean value of the deviations. 

The ‘highest mean’ is entered on the wind record 
simply by visual estimate as the mean value of the 
maximum and minimum speeds over a sampling 
length representing one hour. 

The corresponding record of the deviations from the 
mean speed over the same sampling length forms, 
owing to the slow rotation of the instrument drum 
and the thick trace of the recording pen, a more or less 
irregular jagged area which is determined by the use 
of a planimeter. Division of this area by the scale 
length of one hour then yields an approximate value 
of the mean amplitude of the deviations. 


5. Practical Application 


In order to demonstrate the method evolved for the 
determination of the ‘design wind’ for any particular 
area, the following example will be helpful. 

On the 13th November, 1947, the City of Durban 
experienced a wind storm which caused structural 
damage estimated at well above £100,000. Only once 
in Durban’s previous recorded history, namely in 
May, 1905, did the wind reach similar gale force. 
The anemogram, recorded at the Durban Meteorological 
Office, is shown in Fig. 4. 

It will be seen that the wind rose at about midnight 
and at 2 o’clock was blowing from WSW with an average 
speed of 30 mph and gusts up to 44 mph, increasing 
steadily to an average of 40 mph with gusts up to 
64 mph at 7 o’clock. 

At 7.45 a.m., the wind veered to SW and averaged 
59 mph for nearly half an hour, registering occasional 
gusts up to 78 mph. 





Fig. 4. 














From 8.55 to 9.10 a.m. the gale blew southerly with a 
sustained fury of 62 mph and gusting up to 88 mph ; 
thereafter the wind dropped steadily until evening, 
when increasing heavy falls of rain marked its dying 
violence. 

The period of maximum wind speed and deviations 
over a sampling length of one hour is indicated in the 
enlarged section of the anemogram shown in Fig. 5. 
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Fig. 6. 


It will be noted that the insertion by eye of the 
‘mean line’ presents no difficulty. For the portion 
of the anemogram represented by the sampling length. 
the evaluation yields : 


= 62 mph 
Mean Deviation = 27 mph 


Highest Mean = 
u’ 

Design Wind Uy = (622 + 272)1/2 
G 


_— 


= 67-6 mph 

= 1 + 27/62 
= 1-43 

= 0-0026 G U2 
= 17 lb/ft.? 


The estimated basic wind pressure, ~, refers to the 
design wind, Ui9, computed for the height of exposure 
of the anemometer nozzle, which is normally placed 
at a height of 10m (33ft.) above ground level. 
Depending upon the shape of the structure exposed 
to the wind, the basic wind pressure, ~, must be 
multiplied by the respective shape coefficient, S, in 
order to obtain the wind pressure, P, to be used in 
the design. 


Gustiness Coefficient 


Basic Wind Pressure = p 


6. Variation of Wind with Height 


A factor of great practical importance in structural 
design is the increase of wind pressure with the height 
above the ground. Observations of the flow of liquids 
in pipes have shown that in turbulent flow the velocity 
profile can be represented by a fractional power law 
of the distance from the pipe wall : 

“% = 21/x ° ° ° . (11) 

Experiments which fitted the fractional power law 
were made by Douglas Archibald with anemometers 
suspended from kites up to 1,300 ft. above the ground ; 
they resulted in the formula : 

. Uz/Ui0 = (2/10) -25 : ‘ ‘ - (12) 

In this and the following expressions, the velocities 
Uzand U9 occur at heights z and 10 m above the ground 
level. 
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J. Feld,8 in a discussion on wind pressure, quotes 
two formulae due to two American engineers, C. B, 
Wing and R. Beebee, which are stated below, using the 
notation adopted by the present author : 

Uz/Uio = 1 + 0-003 (z — 10) - : ° > (13) 
and 
Uz/U1o=1+0 -004(z — 10) —0 -000003(z — 10)? - (14) 

On the analogy between wind motion and river flow, 
G. Hellman® expressed the variation of wind speed asa 
logarithmic function of the height above the ground, 
similar to the method employed in hydraulic engineering 
for the calculation of the stream velocities at various 
depths, namely : 

Uz/Ui9 = 0-233 + 0-656 log (z + 4-75) + (15) 

The exponent of the power function which, in the 
author’s view, would adequately represent the velocity 
profile for the purpose of structural design calculations, 
is determined as follows : 

Under extreme wind conditions, the increase of 
wind speed at a height of 100 m from the ground 
above that measured by an anemometer at 10 m from 
the ground is found to be only about 33 per cent. 
Plotted logarithmically as a power function, the slope 
of the straight line from the origin to the 1 -33 intercept 
at the 100 m level, where the ratio z/10 = 10, is given by 

log 10/log1 -33 = 8 
which represents therefore the exponent of the function 

2/10 = (Uz/U30)® - ° ° ° > (16) 
This relationship is more suitably written in the 
transposed form : 

U;/Ur9 = (2/10)9-2 : - (17) 

The various formulae which have been put forward 
to express the variation of wind with height, are 
plotted in Fig. 6. 

The values of the coefficients in Hellman’s formula 
and the powers in the interpolation formulae, are 
influenced largely by the coefficient of roughness of the 
ground surface and even more so by the degree of 
turbulence present in the atmosphere. From the 
surface up to a level beyond the influence of surface 
friction the wind backs and increases in speed at a 
rate depending upon the atmospheric turbulence. 

At about 1,500 ft. level the motion of the air is free 
from the effect of surface friction and is maintained 
by the balance of the forces due to the earth’s rotation 
and to the pressure gradient. The momentum of the 
surface wind is less than that of the upper wind, but 
owing to turbulence there is a tendency towards 
equalisation. 

The structural engineer is mainly concerned with 
wind extremes, when the turbulence level in the 
atmosphere is high and the vertical transport of 
momentum is rapid, so that the wind velocity measured 
at the height of the anemometer nozzle approximates 
more closely to that of the upper wind. 

From the above considerations it is reasonable to 
infer that, where monthly and annual averages of 
upper winds are required, the velocity profiles 
corresponding to the Hellman and Archibald formulae 
will serve. For structural engineering purposes, 
however, the steeper curves plotted in the diagram, 
Fig. 6, should be considered when adopting a height 
factor, H = U;,/Ui0, by which the design wind at 10 m 
must be multiplied in order to obtain the design wind 
for the required level. 


7. Wind-Excited Oscillations 


An engineering structure possess a period of vibration 
which depends upon its mass. The gustiness of the 
wind constitutes a periodic force whose period may 
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ee with the natural period of the structure upon 
which it acts. These vibrations are accompanied by 
stresses set up in the structure. When the period of 
the gust is in resonance with the natural period of the 
structure, the amplitude of the resulting motion is 
much larger than that which occurs when the two 
periods differ. Consequently the stresses induced in 
the structure are increased and may rise beyond the 
elastic limit. The collapse, in 1940, of the suspension 
bridge over the Tacoma Narrows in Washington is one 
of the most spectacular failures of this nature. 

The frequency of vibration due to a wind stream is 
directly proportional to the velocity of the wind. 
Taking the simple case of an overhead transmission 
line, the induced frequency, fj, is given by 

fixkofd - , . (18) 
where 

k = Strouhal number, presently to be discussed, 

v = wind speed, and 

d@ = diameter of the cable. 

The natural frequency, fn, of the overhead line 
may be represented by 


1 T 
f= 5/5 ' : ; - (19) 


= distance between nodes, 
tension, and 
= mass of unit length, 


so that in the event of resonance, when /; = fp, the 
vibrations set up attain to dangerous proportions. 

Even in steady winds, rapid vibrations of small 
amplitude which lead to breaks in cable strands are 
possible due to eddies which form behind obstructions 
such as insulators and tower members, provided the 
wind remains steady long enough for the vibration to 
be set up by resonance. 

Vibrations of transmission lines of large amplitude 
are known as whipping, galloping or dancing. The 
period of these vibrations has been estimated at about 
3 seconds, mainly from the cinematographic records 
of the Hydro-electric Power Commission of Ontario. 

The galloping of transmission lines is generally“ 
associated with glaze storms, and the primary cause 
of these large vibrations may be attributed to a deposit 
of ice on the lee side. The wire, having lost its 
circularity of section, acquires a self-exciting character- 
istic which, drawing energy from the wind, stimulates 
vibrations of a cumulative character. 

Wind-excited vibrations have also been observed on 
tall, self-supporting steel stacks. As in the case of 
overhead high tension cables, a stack exposed to the 
wind stream will give rise to the generation of eddies 
on the leeward side. These eddies develop alternately 
on opposite sides of the stack transverse to the wind 
direction and are carried down-wind, forming a double 
tow of evenly spaced staggered vortices, as shown 
in Fig. 7. 
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Fig. 8. 


A criterion for the stability of a vortex trail has been 
enunciated by Th. von Kdrmdn!° who found that in 
most cases the vortices shed from the cylinder faded 
away rapidly and that only for a definite ratio of the 
distance between the two eddy paths, h, and the 
spacing of the eddies in a single path, /, did the trail 
persist over an appreciable distance : 

hil = 0-283 - . : ° , - (20) 

The generation of the vortices is associated with a 
change in the circulation around the stack which 
gives rise to a lift or lateral force alternating at a 
frequency equal to that of eddy formation, namely 


Lift =Lep 7 D . ° ° ° ie (21) 
where 
D =diameter of the stack 
and 
Ie = lift coefficient, which is assumed to vary 


as a simple harmonic function. 

Experimental values for Le, quoted by R. Ruedy 
vary from 0-61 to 1-05. 

The frequency of eddy formation, f, is given by 
Equ. (18) and may be re-stated in the following form 

f =St.v/D ' : , > (22) 
where 
D = diameter of the stack, 
and 
St = fD/v =a dimensionless ratio, called the 
Strouhal number, after V. Strouhal'?, the first 
investigator of this phenomenon. 

Although the Strouhal number for bluff bodies with 
sharp edges or corners is little affected by atmospheric 
turbulence, this is not so in the case of cylindrical 
objects, as demonstrated by the curve of Fig. 8. 

For atmospheric conditions with Reynolds numbers 
below 10’7, a mean value of St = 0-22 may be used 
in the calculation of the frequency with which the 
vortices are shed. For values of Re approaching 10’, 
atmospheric turbulence reaches a level of intensity 
at which the generation of vortices ceases to be periodic. 

It is clear, therefore, that extreme conditions of 
atmospheric motion and turbulence are less conducive 
to setting up periodic transverse forces which may 
coincide with the natural frequency of the stack. 
This is borne out by the failure of a welded stack of 
16 ft. diameter and 300ft. height, reported by 
Den Hartog8, which came to resonance at its natural 
frequency of about 1 cycle per second in a wind of less 
than 50 miles per hour. 
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Stack vibrations of lower frequency than the above 
have been reported by K. Déring!4, who found that the 
natural frequency of sway of 100m high reinforced 
concrete stacks is of the order of 25 cycles per minute. 
Gust frequencies of this order are not unlikely ; they 
can set up a buffeting action and cause a period of 
resonant oscillation which may lead to the ultimate 
fatigue of the stack materials. 


8. Conclusion 


The complexity of the structure of the natural 
gusty wind is evident from the brief survey presented 
in this paper. There appears unfortunately to be no 
short method of accurately estimating the wind effect 
for engineering purposes. It has been said that “ the 
ideal way would be for the engineer to obtain from the 

eteorologist the maximum wind velocity a proposed 
structure should be able to withstand and how often 
that velocity would be expected to occur.” 

The author’s method of assessing wind extremes 
for engineering purposes merely extends the routine 
evaluation of anemograms at meteorological stations, 
but it is hoped that it may help to discourage the use 
of arbitrary rules in the design of exposed engineering 
structures. Many problems which need further 
investigation suggest themselves from the paper. 

Finally, the author wishes to make acknowledgement 
to the British Electrical and Allied Industries Research 
Association (E.R.A:) for permission to publish this 


paper. 
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Discussion 


The Council would be glad to consider the publication of 
correspondence in connection with the above paper. Communi- 
cations on this subject intended for publication should be 
forwarded to reach the Institution by August 3lst, 1961. 








Book 


Calcul Rapide des Poutres Continues, by M. M. 
Marcel and Andre Reimbert. (Paris : Edition Eyrolles, 
1960). 260 pp. 8}in. x 6in. approx. Price 340 N.F. 
plus postage 35-45 N.F. 

The original theorem of three moments may now 
not be used generally so much in this country, owing 
to the advent of newer alternative methods. However, 
the three moment theorem is fundamental and has 
many applications and uses of its own and cannot be 
dispensed with. On the other hand, the method is 
often long and tedious, so that an approximate, but 
sufficiently accurate, variation to ease the work and 
reduce the time in the solution of cases is always 
welcome. 

This book, compiled and written by two French 
Consulting Engineers, goes thoroughly and deeply 
into cases that arise under uniform and point loads 
on series of varying spans. It can confidently be 
recommended to French-reading engineers for use in 
solving their problems. 

The method follows the rule enunciated by M. Caquot, 
whence, by simple substitutions, the results obtained 
give satisfactory solutions confirmed by tests on 
reinforced concrete beams and slabs. 

The rule is clearly explained in the pages following 
the Preface by M. V. Bollard, who states that it is 
the best and simplest method for these calculations, 
and he compliments the authors on their important 
work of direct application, and for their derivations 
of formulae for spans up to four and more in number 
under different loadings. 


Reviews 












The mathematics are clear and sufficient, while 
the 157 tables, charts and diagrams are adequate 
and give all the essentials required for design. 
Moreover, in worked examples the final application 
to the design of reinforcement is also given. 

LF. 


Reinforced Concrete Piling and Piled Structures, 
2nd Edition, by F. E. Wentworth-Sheilds, W. S. Gray 
and H. W. Evans. (London: Concrete Publications, 
1960). 9in. x Gin. 149 pp. 18s. 

The second edition of this book contains in a concise 
form much information for designers, constructional 
engineers, and students concerned with the manu- 
facture, driving, and design of precast concrete piles, 
the testing of soils, and the design of piled structures. 
The design data is in accordance with the latest 
British codes of practice, and examples of methods 
of design and construction, and illustrations of many 
structures are given. Extensive revisions have been 
made to the previous edition, subjects which are 
additional or have been amplified including loading 
tests, the calculation of bearing capacities of piles, 
the design of sheet-pile retaining walls, and specifications 
and quantities. The stresses occurring during driving 
are dealt with in detail, and methods of designing 
pile-caps for single piles and groups of piles are given 
with numerical examples. The descriptions, working 
drawings, and illustrations of works in the course of 
construction, justify the addition to the title of 
“ Piled structures.” 
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SYNOPSIS 


The behaviour of a structural member subjected to 
a torsional load may be defined by its torsional rigidity 
and the stress distribution over its cross section. As 
distinct from the solid or hollow cross section which 
is homogeneous in the axial direction, the semi-closed 
box section represents a three-dimensioned problem 
of a complex nature. This Paper attempts a solution 
by the Strain Energy method leading to expressions 
of torsional rigidity of a particular type of design and 
also analyses the stress distribution, the results being 
expressed in terms of the applied torque. 

-Considering these factors, a method has been devised 
of designing such girders subjected to combined 
bending and torsion within the elastic range. 

An experimental programme has been carried out 
on full scale girders subjected to pure torsion and 
results were found in general to be in agreement with 
the theoretical work. 





















NOTATION 


cross sectional area of channel 
Ay cross sectional area of bracing rod 
Ao cross sectional area of end plate 

C torsional rigidity 

D° distance between centres of channels 
D, distance between intersections of centre 

lines of bracing rods 

E Young's Modulus 

G Shear Modulus 

H overall width of channel 

I moment of inertia of bracing rod 
moment of inertia of channel about y, 
z axes 

moment of inertia of end plate about 
Rot ype axes 

ength of bay 

force causing torque 

shear force 

moment arm 

applied torque 

strain energy 
shear force 
height of channel 

developed length of channel strip 
distance between bracing rods 
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*This paper is based on part of the thesis submitted by one of the 
authors (H. Nolle) for the degree of Master of Engineering 
Science at the Engineering School of the University of Melbourne. 
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On the Torsional Stiffness of 
Box Section Girders* 


by F. Laszlo,t Dr. Ing. (Darmstadt) and H. Nolle,j+ M.Eng.Sc. (Melb.) 











by width of end plate 
¢ thickness of channel strip 
Cp _ thickness of end plate 

d, diameter of bracing rod 

e = 3(Di — D) 

ly equivalent length of bracing rod 

) angle of inclination of bracing rod to 

axis of girder 

@ angle of twist 

6 angle of twist per unit length 

Poisson’s ratio. 


a 





































Fig. 1. 


1. Introduction 


One of the problems met in the design of steel 
structures has been quite often the determination of 
their elastic torsional rigidity. This has been the 
question not only in cases where members are subjected 
to externally applied twisting moments but also in 
problems connected with lateral buckling in bending. 

Perhaps the first to give a solution to the above 
problem for a framed structure was Batho!, who 
developed a formula for the shear in a structure 
subjected to torsion. For unit loads applied to the 
four corners (Fig. 1), the expression for panel shear 
was given by 


bl 
sa 2 (area ABC) (1) 


The main shortcomings of this theorem were that 
even though the stresses in the members of the lateral 
system are given directly by (1), to find those in the 
main truss members, additional analytical or graphical 
methods had to be applied. Also, the theorem assumed 
the ends of all members to be ball-jointed and subjected 
to tension or compression only. 
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The latter deficiency is not serious when the ratios 
of bending and torsional resistances to length of 
member are small, but in certain cases for more 
compact sections such as considered in this Paper, 
effects due to bending and twisting of members have 
to be included. 

A similar problem, but in a different field, arose 
in the early 1920’s with the advent of the aeroplane. 
It concerned itself in particular with the fuselage 
subjected to torsional loads. Analytical methods of 
solution were developed by several workers?,3, all of 
which, however, involved the second limitation 
mentioned above. 

During more recent years the trend has been 
towards lighter construction of structural frames. 
Among other innovations, efforts have been made to 
replace heavy rolled sections by fabricated sections 
of light weight. The shape of such sections is manifold, 
but a number of these may be grouped under what is 
now accepted as the thin walled fabricated box section. 

The Paper develops solutions to one of the problems 
of stress analysis of such semi-closed box sections, 
namely their torsional stiffness. Solutions are given as 
to its stiffness, strength and load capacity when the 
section is subjected to pure torsion. 

An experimental programme on full scale girders 
showed good agreement with the theory developed. 


2. Method of Analysis 


The type of section under consideration is made up 
of two channels braced on the sides by continuous 
bent rods, the nodes of which are welded to the inside 
flanges of the channels as shown in Fig. 2. The ends 
of the girder are fitted with plates welded to the ends 


of the channels. Such a girder, subjected to twist, 
represents a statically manifold indeterminate struc- 
ture. For low values of torque the strain energy of 
the whole girder may be expressed as a quadratic 
function of the forces and moments involved. Since 
the statically indeterminate forces and moments are 
all internal ones, their correct values must make the 
strain energy a minimum according to the Least 
Energy Theorem. Thus, the partial derivatives of 
the strain energy U of the whole girder with respect 
to every statically indeterminate moment and force 
must become simultaneously zero. This leads to as 
many linear equations as there are statically indeter- 
minate forces and moments, from which the value of 
all statically indeterminate unknowns is obtained. 

The angle of twist of the girder as a whole is obtained 
by Castigliano’s Theorem and hence the torsional 
rigidity. 

A convenient approach is to replace the finite 
girder as shown in Fig. 2 by bays of an infinite girder 
of similar shape and by a finite girder having only one 
bay and with ends as shown in Fig. 2. 
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To find the strain energy of the infinite girder, that 
of one pitch, i.e. that of section AC is to be expressed. 
Though some stiffening is caused by the welds between 
channel and bracing at sections A, B and C, it is 
assumed that over an elementary length Ax at A, B 
and C only the channels transmit the external torque. 

Let the girder be intersected at AA, and the infinite 
section to the right of AA be replaced by forces and 
moments Fxo, Fyo, Fo, Mxo, Myo, and Meo applied 
to each of the cuts across the two channels as in Fig. 3, 
Next the two bracing rods Ja and Jd are separated at 
their junctions with the channels and forces and 
moments Fp, Foyt, Foz, Mv, Mot and Mpz applied 
as shown in Fig. 4. 

The elementary length Ax of the two channels 
take over all duties at BB. For the length BC the 
rods IJa and JJb take their share in the torque. 
Separating the rods from the channels as before, 
the rods are now considered to be subjected to end 
forces and momerts Fy’, F’pt, F’pz, M’p, M’pe and 
M'pz as in Fig. 5. 
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Expressing the strain energy of one pitch and 
differentiating it with respect to the external torque 7, 
we obtain 


oT = Pao 
and 
2LT 
Cac= — 
A Pac 
where 2L denotes the length of one pitch. 
The angle of twist is obtained in a similar fashion 
for the girder having one bay only. Denoting the 
total strain energy of channels, rods and end plates 


PE 
Consider the twisting angle @ of the whole girder of 
length (mn + 1) L due to T, viz. 
_*LT | LT _(n+WLT 
rts Ce C 
from which 
_ (n+ 1)CaCe ; 4 . 3 
c Ps Cac “TT nC py ( ) 


2.1. Equivalent Length of Rod 

For the purpose of analysis an equivalent length /» 
of straight rod replaces the continuous bent rod. 
Two definitions for this length of rod may be considered. 
_(i) Equal Volume—This defines the length of the 
ideal straight rod as such as to make the volumes of 
the bent rod and that of the straight rod equal. 

(ii) Length of Centte Line—The equivalent rod is 
extended to the points where the centre lines of 
adjacent rods intersect. 

In engineering practice the second definition is 
generally accepted; however, solutions by both 
will be given for the purpose of comparison. 
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3. The Girder of Infinite Length 


_It has been shown’ that adjacent bays do have a 
similar stress distribution at corresponding points, 
viz. /Xi3/ = /Xj2! where Xi; and Xj2 are unknowns 
for the first and second bay, respectively. Thus the 
problem reduces to the solution of any bay along the 
length of the girder. 
The equilibrium equations of rods, and channel 
of appropriate periodic length, give 
Foz = 0 ; : 5 ‘ ) 
Feo =0 ; ; ; . , (5) 
Fyo = Fp Sin v Fot cos v . . ) 
Moz = oF vt < “oif2} 


age ) 
Mxo = ul f Feo 2 ‘ . $ (8) 


ae 
LFro = 2Mpyt cos 


2M p sin v 4 
Fytb sin vy - . 


Fybcos 


3.1 Strain Energy 

In the evaluation of half the strain energy per bay 
of length L we have to account for all direct acting 
loads as well as bending and twisting moments and 
shear forces. Localized stresses and stress concentra- 
tions and their effects are neglected, and the channel 
will be considered to be subjected to loads acting on 
its cross section as a whole. 

The bending and tension strain energies are found 
in the conventional way. 

The shear strain energies are for channel# 

: Al+ywLlL,, 
Urro = Ft © (10) 

due to Fyzo. 

Furthermore, the resolutes of rod forces Fp and Fpt 
in the x-direction act on both flanges of the channel 
section and in opposite direction. Replacing the 
concentrated reactions at the panel points by a 
uniformly distributed shear force Q as shown in Fig. 6 
we have 


Q) = ; (Fycos) + Fptsiny) 


It may be shown by means of the membrane 
analogy® that the strain energy of a channel of length 
L is given by 
HQ2L 
a aC 
substituting for Q and G, obtain 

H(1 + p) 


“LE (Fpycosl) + Fytsiny)? 


Ug = 














166 


The only shear strain energy in the bracing rod‘ 
is due to F pt, viz. 


Ury = a i ise 32 , , - (12) 


Denoting the total strain energies of the channel 
and rod by U; and Uz, respectively we have 


L 2 
U; = fe |[~ Fycosy —— Fusing | 





2A,.E 
3(1 + w)L D 
( oUt [F- -% Fro — Mycosy 


2 
— Mysiny — Fog sing + Fus$ cost | 
L 
+ set) [1-— Mypsind + Myrcosy 
°o 


- Fy cosy + Fu: dsing—xF ro] dx 


+a [arse + Foe (3— siny) 


2El, 
2 
9(1 L 

— Fy(e cosy) + ater Fro? 





1-+y)H  #- 
+(Fer [ Focosy + Frsiny | - (13) 
and 
ye 3, i+ (1 + wlpMy? % 
Us 28 + —aer~— e 
Fy(®—« ) fae 4+ 2 uw 
ot\ > — + st vt? 
CT eg = 8 em 
+ 3A-E Fot (14) 


The total strain energy U = U; + U2 as given 
above is that of a girder length of L/2. Remembering 
that Fzo =/(Fp, Fot, Mn, Mont) a solution is given 
by the partial derivatives of U with respect to the 
statically indeterminate quantities equated to zero, 








aU 

aMyo ’ 

dU 

iM. ° 
aU 
~—_ an @ 
OF » 
_ : - (15) 
OF 
dU 

am, ~° 
aU 

OM “ 

The first and second of scat lead to 
Myo=0 ; * (16) 
20 =0 . Z ; - (17) 


4. Torsional Rigidity C 


The torsional rigidity can now be found by the 
principle of conservation of energy of the system as 
a whole. 


For an elastic system, the decrease in 
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potential energy of the external loads must be equal 
to the increase in strain energy due to these loads 
4T OL = UL = 2U 

where @ is the angle of twist per unit length and U, 
is the mean strain energy per unit length of gi: der 
and U is the strain energy over the length L/2 as 
given by (13) and (14). 

Also, by definition 6 = T/C where C is the torsional 


rigidity ; substituting for 6 above, obtain 
fig wea 
E~ 40E a. 


AF 
Pe an 
Figg me | 
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Fig. 7. 


5. The Girder of Finite Length 


The solution of the finite length girder follows the 
same line as that of the infinite girder. Figs. 7 and 8 
show the coordinate system of a one-bay girder with 
solid end plates welded to the channels. Both plates 
have been separated from the channel ends and the 
tentatively assumed direction of torque and internal 
moments and forces between plate and channel are 
indicated. The sense of the forces and moments 
between bracing rods and channel bear the same 
relationship to the sense of the externally applied 
torque as in the infinite girder. 

Applying the conditions of equilibrium, the following 
equations are obtained :— 


Foz = 0 (19) 

Feo = 0 $ « : (20) 
7 if 

xo = 5 —@ Fro é ‘ - (21) 


Fyo = Fy sin Y — Fy cos ‘ : (22) 


l 
Mp: = = Fot ; e ‘ - (23) 
L b y 
Myo = 3 Fro— Fo 5 cos’ — Foe 9 sin 
+ My sinY—Mytcosy - : » (24 


5.1 Strain Energy 


The strain energy expression of one channel is 
slightly modified and reads now 


. L 
ve eee 


— Fycosl) — Fyysinw)2 
Y oerey 


311+ywL/T D b ; 
+ sod Cae Fao + Foe 5 cosy 





TI 
the s 

Th 
of th 
One) 
T : 


expel 


Thus 


Now 
girder j. 
(27) an 





M pcosy — Mysin ») 


» Boe 
Fy5 siny 


2 
x) Fuzo2dx + —_ 
ail; 


e sin +) — Fyecos | 


4EcL 


9 
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(Fycosy 
(25) 
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Fig. 8. 


The strain energy expression of the rod is again 
the same as before, 1.e. equation (14). 

The total strain energy of the plate will be a function 
of the external torque 7 and its mode of application. 
One way of applying a torque is shown in Fig. 8 where 
T= P.R, and conforms with the loading of the 
experimental specimens. 

Using the same method as before, the shear strain 
energy in the plate due to shear force Q is given by 

- 1+ yn)AB _ ; 
Ug = TEE (Fp cos Y + Foyt sin )? 
(26) 


Thus the strain energy Us of the plate becomes 


. D ae : 2 
U3 = 9AE ( psinu F wcosi) 


, A +p)D f(T F z 
8A ,E R 7 


3(1 + n)D/L_ > 6 
ars ( Fao — Fg cos 


« 
~ 


«ew 
—Fore5 sin 


. D 
+ Mysind — Myxcos ») + oT 
& min 


D fe 
(2-1) 


Mic? 


(l+ywAHB_ : wee 
+ SE cE - (Fycosy + Fryesiny)? 
_Now the total strain energy U of length L/2 of the 
girder is given by the sum of equation (14) (25) and 
(27) and the simultaneous solution of the partial 
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derivatives of U with respect to the statically indeter- 
minate quantities equated to zero leads to the remaining 
unknowns. Thus 


aU 
OF p 
aU 
OF vt 
aU 
dMy 
oU 
OM pt 
oU 
@Mzo 
oU 


5) aa 


The torsional rigidity is given by (18) after evaluating 





U. 

It may be noted that the solution of the latter type 
i.e. one for a single bay girder may be used to compare 
experimental and theoretical results directly without 
the help of equation (3), i.e. it eliminates uncertainties 
arising out of superposition. 

In Table 1 solutions of equations (15) and (28) are 
shown for three girders. The dimensions used in the 
calculations correspond to those of the specimens. 
Torsional rigidities have been calculated with the two 
different definitions of equivalent rod length as 
mentioned before. 


6. Structure Assumed Pin-jointed 


The foregoing method involves a considerable amount 
of numerical work in the solution of equations (15) and 
(28). A largely simplified, but less accurate answer 
is given by treating the bracing as tension or compression 
members simply supported at both ends. End 
moments and transverse forces become zero and Fp 


acts alone, viz. 
My = Mut = For = 0 
Note that the ball-jointed structure in this case 
represents a similar problem as dealt with by Batho 
(loc. cit.), however, the structure is still statically 
indeterminate and requires the Least Strain Energy 


method. 
With equations (9) (13) and (14) the statically 


indeterminate quantities Myo and Fp are obtained 
from 
aM,, = ° | 
yo ° : ; a 
aU _ (29) 
OFp 

It may be shown that My, is zero again. 

Similarly to the infinite girder, the finite case is 
solved by the partial derivatives of U with respect 
to the unknown quantities 

aU _ 

OF yp 

OMzo 

<a 

OF zo 
where U is given by equations (25) (14) and (27). 

Values of torsional rigidity by this method are 
shown in Table 1 for three girders. 
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TABLE 1. 
Theoretical values of statically indeterminate quantities and torsional rigidity for girders A6, B6 and Al. 





Girder A6 B6 





Type of 
solution | Structure rigid Structure Pin Structure rigid Structure Pin Rigid Pi 
jointed jointed joint 








Equation | (28) (29) (30) ( (28) (29) (30) (28) (31 





Fy/T 0-15511 0-15129 0- 15689 0+ 15450 0-14779 0-14452 0-14916 014623 0-15456 | 0-15 








Fp,/T 0-00017 0-00019 0 -00002 0-00004 0-00020 | 





My/T |—0-00109 |—0-00129 —0-00038 |—0-00094 0-00135 | 


Myt/T 0-00307 0 -00333 0-00275 0 -00307 0-00355 | 


Mzo/T 0-00073 0-00076 0-00081 0-00016 “00071 | 0-0 





| 
Fzo/T 0-02661 0-02625 0 -02347 | 0 -02837 -02663 | 
| 


0-02 351 





3-992 4-029 5 | 3-840 
3-763 3-786 | 3-61! 


4-79 3-798 
C/E 4-59 3-573 











CoolE 





With ‘ 
eq. (3)Cp/E * 325 : of 5-247 
. . . +997 











Experimental 
C/E 














Note : Upper figures for C are calculated using definition (i), lower figures using definition (ii) of equivalent length 
of rod (see section 2. 1.). 


Te 


7. Non-Linear Characteristics UE the 
can not be expressed in the form 7 is not pro- sectic 


The experimental curves for torsional rigidity . en 
which were obtained from each girder show clearly portional to T?. : 
that C does not remain constant even in the range Instead of U we may write (U’ + Uc) where U’ Th 
where all internal stresses are below the yield point. is now the total strain energy due all forces and the | 
This is due to the non-linear behaviour of the bracing moments with the exception of M, having the form mach 
rods subjected to simultaneous axial compression kT? where k is a constant, and U¢ is the variable strain a 
and bending. These are the rods under compression energy due to the non linear M. The latter can also be - ; 
and bending due to Fy and Mpz and Mpt, respectively. U 1S. : 
For practical purposes, the effect of Mpz may be reduced as =5 for any value of T. peay 
neglected since Myz < Mpt and Mp; acts to produce T specit 
double curvature whereas M pt causes single curvature. The modified expression for torsional rigidity is with : 

The bending moment at any point x from the end given by way 1 
of the rod may be expressed by an approximate P LT?2 of the 
method® frequently used in the solution of problems E~ 40° + UE . > ( Cur 
involving combined axial and lateral loads. Thus , ‘ =a torqu 
the bending moment is This can be expanded as a power series in 7 or Stre 

calculated for individual values of T, the latter being leas 

V= Melt F y(lx — x*) } P preferable because of slow convergence of the series _ 

J - M pt a i FIL (31) | 

2EI(1 a) at large values of T. jemnr 

Fyl2 Substituting equation (32) into (3), and remembering Stress, 

where a = aan and / is the length of the rod between that only half the rods are under compression, equation ae 0 
n?EI (3) becomes I 

ne 


pin-ends. For the purpose of evaluating M, the Cc strain 


length is not equal to our previous /p». It will be ee henlic 
assumed that the ends of the rods acting as ee (n + 1)LT2 a 
are located at the mid-width of the channel flanges me 7 Tia) A 
for the infinite girders and for the single-bay girders / 4E[n(U% —4U® yott+3Uc® ) + (Ui—}U mot + 31 ea) term 
terminates at the edges of the channel flanges since . ; 33) Th 
in the latter case flexure of the flange is prevented by where Ujypt and U® yp¢ are the strain energies due mom: 1 
the end plates. to Mp, for the finite and infinite case, respectively. torque 
The bending moment is no longer proportional to C—curves, based on equation (33) have been drawn values 
the applied torque and consequently the strain energy in Appendix 1. tods. 








8. Experimental Investigation 


Tests on full size girders were carried out to study 
the torsional properties and stress distribution in box 
section girders discussed in the Paper. The dimensions 
of each girder are shown in Appendix 2. 

The size and length of the girders were limited by 
the load capacity and length of lead of the torsion 
machine, respectively. As the maximum length of the 
specimen that could be accommodated was six feet, 
two nominal lengths of six and one foot were chosen. 
Fig. 9 shows a 6 ft. girder fitted into the machine and 
ready for testing. A pure torque is applied to the 
specimen through a pair of link mechanisms provided 
with spherical roller bearings at all link joints. In this 
way it was possible to apply a pure torque irrespective 
of the deformations of the girder. 

Curves of torsional rigidity plotted against applied 
torque are shown in Appendix 1. 


Stresses were measured by strain gauges attached 
to various parts of the girder as shown on Figs. 10 and 
ll. Points of interest were mainly bracing rods under 
compression and, to a lesser extent, rods under tension. 
Stresses in the channel section were measured at 
*=(0 and x* = L/2. 

In‘ernal forces and moments were calculated from 
strain readings corresponding to an increment of the 
applied torque of 40,000 in.-Ib. 


A complete set of forces and moments expressed in 
terms of T is shown in Table 2. 


The non-linear changes of the internal forces and 
momcnts when measured against increasing applied 
torque is demonstrated in Figs. 12 and 13 showing 
of Fp/T and M»,/T for tensioned and compressed 
rods. 





Fig. 10. 


8.1 Mode of Failure of Girders 


Girders A6 (6 foot) and 1A (1 foot) were tested to 
failure whereas girder B6 (6 foot) was loaded only 
within the elastic range. 

(a) Test on Girdey A6—Failure of the girder was 
initiated by simultaneous buckling and yielding of 
rods under compression. The first rod to buckle was 
one connected to the centre of the girder and yielding 
commenced at 75,400 in.-lb. 

After some further twisting, a compression rod 
midway between centre and end of the girder buckled, 
followed immediately afterwards by a rod reaching 
the end of the girder. The ultimate load recorded was 
79,100 in.-lb. The deformed girder after a twist of 
approximately 18° is shown in Fig. 10. 





The Structural Engineer 


(b) Test on Girder Al—First signs of yielding were 
observed at 85,000 in.-Ib. The yield was due to direct 
overstraining and no signs of flexure or buckling could 
be seen. A considerable amount of plastic stain 
developed at a steady rate at approximately constant 
load before inwards buckling of the bar took plice. 
The ultimate load reached was 95,100 in.-lb. 

The specimen after the test is shown in Fig. 11. 


8.2 Experiment Versus Analysis 


An examination of Table 2 shows that a consider: ble 
variation exists in the magnitudes of the forces ind 
moments measured on different rods and parts of the 
channels. Two causes have been considered to be 
responsible. 

(i) Lack of uniformity of the girder itself—All dimens ons 
in Appendix 2 are averages. In fact, appreciable 
deviations were observed in the angle nd 
dimension 6 which is most important as far as 
stress distribution is concerned. [Initial curvaiure 
of the rod, either inwards or outwards, will tend 
to increase or decrease respectively, the benc ing 
moment in the rod under compression. 


Fig. 11. 


TABLE 2. 
Experimental values calculated from Strain Gauge readings. 





Mp Mot 





GIRDER A6é- | Fy Foe 

| Rod = =—Ss—~CS~*=éd a | (01332 *| 000057 —0-00141T | 0-00928T 
014287 

“Illa | 0-13887*| 0-000467 0-00058T | 0-0148T 
| O- 


0-00058T |} 0-0122T 


Va | 0-136 


IlIb | 0-1300T 0-00575T | 0-0128T 
Ila | 0-12697*)} 0-00066T —0-00406T | 0-0060T 
-1438T 


0-0115T 





Channel section | | 0-07747¢| 0-101T+¢ | 0-0137 
| . ol me a nn = > ——-——} H 
3—3 | 0-0042T 0-100T+ 0-937T+ | 0-0057 


| (O)x - 

















4 
GIRDER Al 
: ae) Ee Ee. A 
Ia | 0-1253T*| 0-00068T 
0-1392T 


Ib | 0-1233T*| 0-00069T —0-00645T| 0-00995T 
0-1590T 
Channel section 1—1 0-099T+ 


GIRDER B6 


Siti WS 
| 


—0-00407T| 0-0108T | 


























“ta! 0-127T | 0-00022T 0-0132T 





_ - Ea a 





Illa | 0-126T | 0-00038T 0-0176T 
aed ee ae si 
0-C178T 





IIb | 0-127T | 0-00030T 


‘ten “IVa | 0-129T 0-0121T 
Channel section 3—3 0-1187 —\-flange gau 
0-141T —/-web _,, 














PR I I's af 0-0827 —\|-web 









































* refers to strain gauge readings at the ends of the rods, unmarked values are from gauges at midlength. theref 

er 
+ these values may be + depending on whether Fp is tensile or compressive. we 
All values correspond to 7; = 5000, Tz = 45000, A T = 40000 in. Ib. bracir 


Note. 
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(ii) Localized and non-uniformly distributed stresses—In 
calculating the internal forces and moments the 
assumption was made that all forces and moments 
were acting through centres of gravity and about 
neutral axes, respectively. This is not valid in 
this case. All forces and moments applied by the 
rod to the channel act through the weld at the inner 
urface of the flanges; thus axial direct stresses 
will not be uniformly distributed and this is 
demonstrated by the values of Fp and Fy, for 
girders A6 and B6 in Table 2 (Fx = — Fpcosy — 


Fyesin ). 
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Fig. 13. 


The deviation from simple beam theory type bending 
stress due to moments applied by the rods to the channel 
is possibly of a similar, if not larger magnitude. The 
presence of high local bending stresses was established 
beyond doubt by additional strain measurements. 

In conclusion it may be said that the measured 
values in Table 2 which refer to the channel cannot 
therefore be compared with the analytical results. 
However, a comparison is more satisfactory for the 
bracing rods. 














8.3 Non-Linear Torsional Rigidity 


The graph showing C/E versus T in Appendix 1 
indicates a gradual decrease in torsional rigidity until 
yielding or buckling load is reached at which stage the 
collapse of the girder is marked by a sharp increase in 
slope. 

The assumption of constant ratios versus torque T 
of the internal forces and moments (with the exception 
of M») is confirmed by Fig. 12 which shows a variation 
in F»/T of only 6 per cent over the range of T. 

The considerable variation in Mp¢,/T on the other 
hand, is seen in Fig. 13 where mean values of M»p,/T7 
are plotted for compression and tension rod. 

An interesting feature here is the fact that a larger 
rigidity (C); = o does not necessarily imply that the 
reduction in C with increasing JT will be at a smaller 
rate, nor does a greater initial rigidity mean that a 
higher ‘ ultimate load ’ will be reached before buckling 
took place. 


9. Application of Analysis to Design of Girders 


A girder such as considered in this Paper, rarely is 
subjected to torsional load alone. As a structural 
member, the normal loading conditions are most 
likely to be such as to cause a bending moment and 
possibly axial compression. 

From simple beam theory the bending moment 
and shear force capacities are given appropriately by 

M = cywA,D 
V= (Ft — F,) sin v 
where cw = working stress in the compression member 
D = distance between centres of gravity of 
channels 
F;, Fe = maximum permissible forces in the 
tension and compression rods, 
respectively. 

The length of the bay is controlled on the one hand 
by bending requirements. For least weight, the bay 
may be made of such a maximum length for which 
the compression ‘flange’ is at the limit of elastic 
stability. This length and D will define the shear 
capacity of the girder for any rod size. The sections 
treated in this Paper have | = 60° which corresponds 
to a design for bending moment and at the same time 
maximum torsional rigidity (see Appendix 3). 

When shear strength is of primary importance, 
maximum rod efficiency will be attained when the 
rods are sloping at 45°. 


9.1 Eccentric Load 

A twisting moment will be set up when the line of 
action of a transverse load does not pass through the 
centre of gravity of the section. Twisting of the section 
will take place until equilibrium is reached. The 
torque is now the sur due to initial eccentricity and 
eccentricity caused by twisting of the girder. The 
twist is calculated using the C value corresponding 
to equilibrium torque. 
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The true stresses in the girder under combined 
bending and torsion can be found by superposition. 
Let the shear force be V, corresponding to a bending 
moment M, then at section a-a, Fig. 14, the shear is 
resisted by two bracing rods. Provided the girder is 
not overloaded, F,t = Fe = F and V = 2F sin | or 

V 

~ 2Qsind. 

The force in the channel at a-a is M/D causing 
tensile or compressive stress. 

If a torque T is superimposed as shown at a-a, a 
compressive force «aT = F’ is introduced in both rods, 
and the resultant force in the rear rod is (F, + F’) 
and for the front rod (F;— F’). In the next bay the 
rear rod is under tension of (Fy + F’). Thus a non- 
symmetrical stress distribution is set up in the bracing 
with maximum and minimum values in the rear plane 
of rods. 

Collapse of the girder will mitiate in the rod under 
compression of (F, + F’) or, substituting from above 


(ear + 27) where « is the constant obtained 
2 sin 


from the analysis of the girder subjected to torsion. 

In the particular case of large bending moment 
and torque, the combined stresses in the channel 
ought to be checked as well as in the rods. The axial 
force in the channel due to 7 changes sign at each 
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panel point and may be additive or subtractive to the 
uniform force caused by bending. 

When checking the strength of a girder, the forces 
and moments obtained for the infinite length girder 
should be used. The difference between (Fp) © «nd 
(Fp), was found to be small, nevertheless (Fp) © is 
slightly larger thus representing the critical case. 

This feature was confirmed experimentally by the 
mode of failure of girder A6 where rods under 
compression and furthest from the ends buckled first, 


10. Conclusions 


. It has been shown that strain energy methods can 
usefully be employed in the calculation of the 
torsional stiffness and stresses in a box section 
girder. 

. The torsional! stiffness of the type of girder investi- 
gated is is not constant but varies slightly with the 
applied load, the maximum value occuring at zero 
load. However, the initial stiffness is not necessarily 
a criterion of its ultimate strength in torsion. 

End effects owing to the method of end finish may 
considerably influence the behaviour of very short 
girders. 

. The critical stresses were found to be in the bracing 
rods, and the structural behaviour under torsional 
load is largely determined by their arrangement. 

. The method, as outlined, in the paper, can be used 
for the design of a girder under torsion and bending. 
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APPENDIX 3 


Effect of Length of Pitch on Torsional Rigidtty 

During this work it was realized that one of 
the major factors affecting the behaviour of the girders 
was the inclination of the bracing rods to the axis of 
the girder. 

Considering the solution as given in Section 6, the 
tod force Fp can be expressed in terms of } by sub- 
stituting in (29) for L and /p 

L = D, cot 

Dy 
sin 
where D, = vertical distance between intersections of 
centre lines of bracing rods. 

Hence, after simplifying, 
Tb, sin Y cos 
b= 3,1 : ch ; v * = (35) 
aycos®) + agcosysin®) + ag 
3(D + Dy,)b(1 + p) 


2ac8 
1 b2 e2 
a4 = ds| 4: + i27, +f] 
D2 
3{ 2D + Di + D; b2(1 + p) 


aa = 
2ac3 


and ly = 








where 6; = 





12(1f+ u)b?  (l+y)H 
5HcD; " 2D; 
_ Di 
ay 8 
and e = 3(D, — D). 
Substituting (35) into the expression for strain 
energy U, we find 
Co 
E 





+ 


43 


aycos8 +agcos psin® )+azg 


dD, 
2 | 2cy(a;cos? b+ agcos sin? )+agz)—b,2cos sin?) 





(36) 
3(1 + p)D; 
4ac 

Curves of Co/E and F»/T calculated fot girder A6 
are shown in Fig. 15. It is interesting to note that the 
values of C% is greatly dependent on in the range 
used in engineering, reached a maximum at approx- 
imately 59°. 

Force Fp, on the other hand, remains fairly constant 
at large values of | but increases rapidly when pp falls 
below 45°. It may be seen that maximum stress in 
the bracing of the girder is reduced considerably by 
choosing a large value of v. 


where c, = 


40 ft 


A | 
1.277 (|3-05°) 


——+ 





4.057 (58-6°) 





Fig. 15. 


ACKNOWLEDGEMENTS 
The work described in this paper was carried out 
in the Civil Engineering Department of the University 
of Melbourne, and the authors would like to thank 
Professor A. J. Francis for making facilities available. 
The authors are indebted to Standard Steel Pty. 
Ltd., Melbourne, Australia for the girders tested. 


References 
1. Batho, C., Engineering, Oct. 15, 1915. 
2. Pippard, A. J.S., and Douglas, W. D., R.&M. 736, Jan. 1921, 
A.C.A., v.1, 1920/21. 
. Southwell, R. V., R.6M. No. 737, A.C.A. v.1, 1920/21. 
. Timoshenko, S., ‘‘ Strength of Materials,” 3rd Ed., v.1, p. 319. 
5. Timoshenko, S. and Goodier, J. N., ‘‘ Theory of Elasticity,” 
2nd Ed., p. 298. 
3. Timoshenko, S., ‘‘ Strength of Materials,’ 3rd Ed., v.2, p. 52. 
7. Nolle, H., ‘“‘ On the torsional stiffness of box-section girders,”’ 
Thesis presented M.Eng. Sc., Melbourne, 1958. 





The Structural Engineer 


A Method of Minimum Weight Design 
for Thin-Walled Beams 


by S. Krishnan and K. V. Shetty 
National Physical Laboratory of India 


Notation 


local buckling stress, p.s.i. 
bending stress, p.s.i. 
yield stress, p.s.i. 
Young’s modulus, p.s.i. 
flange width parameter, in. 
semi-web depth, in. 
ratio h/b 
wall thickness, in. 
stiffener spacing, in. si 
number of stiffeners along each width. 2 2 © Me) 
area of material in stiffener section, in®. 1+ avin'r2 +4) Ls 2bei+K4s) 1 = 40+ KY 
area of material in the beam section in?. Fig. 1a.—Unstiffened Sections 
moment of inertia, in‘. shear! : 
constants for local buckling stress. oo 
Aa, Ag Lagrange multipliers. ae 
M external moment, in.-Ib. 
M,, Mz, Mz reduced moments, in.-lb. 

@ slope change per unit length (M/EJ); 


a 
4 


A= 2b 2+k) A * 2bt(i+ x) A = 4bt(i+k) 


Amma 2 aawrerotnd 2 


Ee 


= HY) 


Introduction 


Design of any structure for minimum weight is Fig. 1b.—Lipped Sections. 
essentially based on the determination of optimum 


proportions on the basis of simultaneous failure in ~~) 
all the critical modes. The critical modes for a thin- “4 o 
walled beam subjected to bending loads would be | 





due to yielding of the extreme fibre in tension or 
compression and local buckling of the compression 
flange. For simultaneous failure 

Cer = Cp =ay - ° . . : - @& 
If only two quantities are to be determined for a 
particular type of beam, equation (1) would suffice as2p[uzeu)+2q4] » *2p[turndg| a osonean*9 
to fix fully optimum dimensions. Apart from the modes 
of failure referred to, special restrictions such as a 1 + PE [ri2+ yee] 1 «2B luirnsr+og] 12 48K [tek )+ 04 
limit on the maximum slope change may have to be ! t 4 ; 
allowed for in certain cases. This means there would see F 
be one more restriction and the conditions for optimum oe Cerne Sacttene. 
design may be expressed as 


Ser = Sb = Fo = Oy Unstiffened Sections : 

and @=4, - . : ; ; ‘a Figs. la and 1b describe unstiffened I-, channel, 
where o@ represents the stress corresponding to and box sections with dimensions and sectional 
the limiting slope change per unit length and @, the properties as shown. In deriving expressions for the 
limiting slope change per unit length corresponding optimum dimensions the following assumptions are 
to the yield stress ay. made :— 

It would also be possible to use these equations to 1. The flange and the web are of the same thickness. 

determine optimum proportions in terms of a structural 2. The web-depth is related to the flange-width by 
loading coefficient! or a structural index?. However, a ratio k which is treated as a variable. 
a direct approach to the problem of designing for 3. The lips, in the case of lipped sections, change only 
least weight may be made by minimising weight or the end condition of the outstanding flanges and 
area of material in the cross-section subject to the do not appreciably contribute to the area and 
restraining conditions expressed by equations 1 and 2. moment of inertia of the section. 

In this paper formulas have been derived for optimum . The flange width 6 is determined by the know 
dimensions of thin-walled I-, channel, and box sections value of the yield stress thus having only one 
with and without lips, and with and without stiffeners. restraining condition and two variables ¢ and 4. 
Curves from which optimum dimensions may be Substituting corresponding terms for ocr and % 
directly read have also been plotted for steel and in equation 1, the restraining conditions for the 
Al-alloy. three types of sections may be written as 
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TABLE 1 





I-Section 


Channel Section 


Box Section 





k=1-3 = 0-646 
t = 0-541[M/KE]}1/3 


b = 0-541(KE)16M1/3g,-1/2 





t = 0-86[M/KE}}/8 
b = 0-86(KE)1/8M1/8c,-1/2 


k = 0-646 
t= 1-084(M/K’ E}1/8 
b = 0-542(K’ E)1/6M1/8e,y-1/2 








2 M 

~ 2ty2k(2 + F/3) 

* M 

~ Dfy2h(1 + B/3) 
M 


en, ‘1 PP —- omnes 
KE (t/2v) tp2h(1 + R/3) 


KE(t/p)? (I-Section) 


KE(t/p)? (Channel Section) 


(Box Section) 
(3) 


While minimising the area subject to the restraining 
conditions (3), Lagrange multiplier 4, is used leading 
to the variational equation 
3(A + A1¢1) = 0 
Where A is the area of the material in the cross- 
section and 9 the restraining condition. The restraining 
condition 9 is expressed as 
91 = Scr — Oy = O 
Differentiating with respect to ¢ and k 


0A Og: _ 
in =” 


aA 001 
2" s 
; ? - (4) 
Substituting for the derivatives of A and 9) in 
equations (4) 
2b(2 + k) + A13ék(2 + k/3) = O 
2bt + Axf3(2 + 2k/3) = 0 
(I-Section) 
2b(1 + &) + Ay3#A(1 + 
2bt + Auf3(1 + 22/3) = 0 
(Channel Section) 


/3) =0 


k 
4b(1 + &) + 243R(1 + 2/3) = 0 


4bt + 24f3(1 + 2k/3) = 0 
(Box Section) 
j — - 6) 
By eliminating 4; from each pair of equations 
solutions for ¢ and & are obtained from the resulting 


020r 


= 


4 
z£ - 
——~ Box on? * 
“Bi steel 
eo 


4 lie —— —$—$_$_—4 


20 30 42 50 








MX 10°Lb-in ——— 


Fig. 3.—t vs M curves for unstiffened Al-Alloy and 
Steel Sections. 


equation together with the corresponding restraining 
condition. The results are shown in Table I. 

Coefficient K for unlipped sections and K+ for lipped 
and box sections are read from standard curves, and 
found to be 0-45 and 3-6 respectively’. Values of 
# and 5 have been computed for a range of moments 
1,000 to 50,000 in.-Ib. for unlipped steel and Al-alloy 
sections and plotted in Figs. 3 and 4. The compressive 
failure stress has been assumed to be 44,000 p.s.i. for 
Al-alloy and 40,000 p.s.i. for steel. Young’s modulus 
E is taken as 10-3 x 10%p.s.i. for Al-alloy and 
30 x 10% p.s.i. for steel. For lipped I- and channel 
sections, the same solutions as for unlipped sections 
would apply except that K will have to be replaced 
by Ki}. 


Stiffened Sections : 

Fig. 2 shows the stiffened I-, channel and box 
sections with dimensions and sectional properties 
marked. The assumptions used in this case are :-— 

1. The web and flange are of the same thickness. 

2. The web-depth is related to the flange-width by 
a ratio k. 

3. The web stiffeners do not contribute appreciably 
to the moment of inertia and they are ignored for 
this purpose. 

. There is a limit on the slope change which is reached 
when simultaneous failure occurs due to bending 
stress and local buckling. 


. The stiffeners are of equal area and uniformly spaced. 


2 | 


Channei Section 





—$——_—_— 4 —E 


30 40 50 


-3 
MX 10 Lb- in ——— 


Fig. 4.—b vs M curves for unstiffened Al-alloy and 
steel sections, 
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Substituting the appropriate terms for oer and op 
in equation 2 and for a limiting value of slope change 
per unit length the restraining conditions for the three 


types of sections may be expressed as 
M 
KOE (t/4)" = SaRe + BS) + Dad 
ans M 
~ E ~ 2y3k2[t(2 + 2/3) + 2a/d) 








(I-Section) 
M 
2p*h[e(I + &/3) + afd] © 
M 
° = 5-5 H + A/3) + ald] 
(Channel Section) 





K1E(t/d)? = 





M 
4p -R((1 + &]3) + a/d] 
M 





K1E(t/d)? = 





° = Fae + AS) + afd] 


(Box Section) 
(6) 


With the aid of Lagrange multipliers Ay and Ag the 
area A is minimised subject to the restraining conditions 
(6) giving the variational equation 

3(A + Argi + Avge) = 0 
Where 9 and gg are the restraining conditions. The 
restraining conditions g; and 92 are expressed as 

91 = Scr — op = O 

os = El —_ 6 = ~- 0 
Differentiating with respect to ¢, d, and 6 

oA 091 oo2 

Ste tee. 

0A o91 O92 

"2°" 


aA 2 


ae ea one _ 9 
b 


+2 3 
(7) 


Substituting for the derivatives of ¢, d, and } in 
equations 7. 


2 
2b(2-+8) +0 Fe BR +R/8) + “7 | te +418) = —0 
—4ab/d2- narra  (2+8/8)-+6a/a* | Agb2h2 -2a/d2?=0 


2[t(2 +k) + 2a/d) + a ae - 2b [t(2 + k/3)+ 2a/d] 


+ A23b2k2[¢ (2 + k/3) + 2a/d]) = 0 
(I-Section) 


TABLE II 
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2b(1 + &) + a - [30 + k/3) + 7 
+ re + b3R2(1 + 2/3) =0 
+ h/3) + sa/ae| 
—rgb3k2a/d? = 0 
- 2b4(1 + k/3) + ald) 


+ k/3) + a/d} =0 
(Channel Section) 


2t 
— 2ab/d® — ry - £202 [ al 


2k 
2[e(1 + A) + ad] + - 


+ dg + 3b%R? [¢(1 


b2k _ 2at 
4b(1 +k) +01 | sec + ms) + a] 


+ Agb3k2(1 + 2/3) =0 
— 4b(1 + R) a di + #252 | a + k/3) 
+ 3a/ds | — reb3k2a/d? = 0 


A(l +4) (t+ a/@) +4 Ss - 2b[4(1 + B/3) 
+ a/d) + dre - 3b7k2[¢(1 + 2/3) + a/d) =0 
(Box Section) 


(8) 


A, and Ag are eliminated from each set of three equations, 
and the resulting equation together with the corres- 
ponding two restraining conditions yield solutions for 
6, t, and d. The flange width } divided by the stiffener 
spacing @ gives the number of flange stiffeners . The 
results are given in Table II. 


The various functions f, g and A of & are as 
defined in the appendix. The results in Table II may 
be grouped together in the following form into one 
set of solutions which can be readily applied to any 
section :— 

b = [K’E — 4a— 26 — 5M, 41/18 
t = [K’ — 8E — 146 §2M)!/18 
n = (K’-2E-5q- 99 - 3M45)1/18 


(9) 


(Channel Section) 
(Box Section) 


Where 
M, = 1: pod lag airy . (I-Section) 
= 1-506 [g:(k)}}/4M (Channel Section) 
= 0-749 tia ) 4M (Box Section) 
Mz = 1-507f2(k)M ; (I-Section) 
= 1- See (Channel Section) 
= 48he(k)M- (Box Section) 
Mz = 1-504[fs(R) }1/5M (I-Section) 


0-7 
1- 

J fn [gs(k) ]/5M 
0-749 [hg(k) }1/5M 





I-Section 





Particulars 





Channel Section 





Box Section 





| 1-134[K’ E-4a-2 6-5g,(k) M4}3/18 





b | 1-134 K’ E-4a-2 6-5f, (k)M4}1/18 0 -915[K’ E-4a-2 6-5h;(k) M4)}/18 
t 1 -032[K’-3 E-1a6 62f2(k)M}1/18 1 -032[K’-8E-1a6 §2g2(k)M)}*/13 | 0-978[K’-3E-1a% 62h2(k)M)}/i8 














n | 1-170[K’-2E-Sa-® 0-3f4(k)M5)1/13| 1 -170[K’-2E-Sa-9 6-3g9(k)M3}1/13| 0 -895[K’-2E-Sa-96-9hg(h) M94 
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Fig. 11.—Variation of M, with M for different types of 
sections. 
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Fig. 12.—Variation of Mz with M for different types of 
sections. 


Values of 6b, ¢ and m are computed for a range of 
moments 10,000 to 1,000,000 in.-Ib. for both steel and 
Al-alloy and plotted in Figs. 5, 6, 7, 8, 9 and 10. 
A value of 0-2in?. has been used for stiffener area 
and the coefficient K’ has been assumed to be 3-6. 
Corresponding values of My, Mz and Mz for a given 
M should be read from Figs. 11, 12 and 13 and used 
for reading b, ¢ and m from the respective curves. 
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Fig. 13.—Variation of Ms with M for different types of 
sections. 


Discussion 


It may be observed from Figs. 3 and 4 that for 
a given bending moment a larger thickness is required 
for unstiffened aluminium alloy section over steel but, 
at the same time, the width could be smaller. Choosing, 
for example, a channel subjected to a bending moment 
of 30,000 in.-Ib. the optimum dimensions for steel 
and aluminium alloy are found to be 

t b area 
Steel 0-110in. 2-04in. 0-742 sq/in. 
Al-alloy 0-158in. 1-62in. 0-846 sq/in. 

It can be seen that for a given moment the optimum 
Al-alloy channel can be over 50 per cent lighter than 
the steel counterpart. 

Again from a comparison of the lipped and unlipped 
sections, it can be seen that K’ being eight times K, 
optimum #¢ for lipped channel is reduced to half that 
for unlipped while the value of b, for lipped channel 
is increased to 4/2 times the corresponding value 
for the unlipped. There is, therefore, an overall 
reduction in weight amounting to nearly 30 per cent 
if lipped channels are employed. This agrees with 
results obtained for channel sections in another 
paper‘. 

It can also-be shown that for the same bending 
moment the optimum J section derived by the above 
method is lighter than that obtained by the structural 
Index Technique described in Ref. (2) °. 

For stiffened sections a study of expressions (9) 
and the corresponding curves reveal that while the 
flange width and number of stiffeners are considerably 
influenced by the moment and the limiting slope 
change, the thickness is not so sensitive to changes in 
bending moment. While using the results, fractional 
figures obtained for shall have to be rounded off to 
and the next higher full number. 


Conclusions 


It should be theoretically possible to solve for any 
number of unknown dimensions by introducing them 
as variables in the function to be minimised and 
imposing more restraining conditions. But too many 
restraining conditions ° may reduce the -practical 
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applications of this method and may not yield reason- of the compression flange takes place only when the 
able results for design purposes. At least one of the bending stress develops the value cy. The method 
leading parameters should be left out of minimisation can be readily extended to the plastic range by taking 
and determined by other physical considerations. suitable functions for oer and op. 

One feature or advantage of this. method is its 

enerality. It has not been restricted to a definite 
ae for cp in the general treatment. This enables Acknowledgement 
the designer to have full freedom to choose his other The work described in this paper is a part of the 
parameters depending on what design stress value research programme of the National Physical Labo- 
he takes. In the present paper the authors preferred ratory of India and is published with the permission 
to select the value for 6 in such a way that the buckling of the Director. 


Appendix 
I-Section 
The equations (8) for the I-Section may be rewritten 


as 
2d2(2 + k) + rrbke[3¢(2 + /3) + 4a/d] 
+ Agd2b2R2(2 + k/3) = 0 
Qad + r,bkt2[t(2 + k/3) + 3a/d] + roadb®h2 = 0 
2d2[t(2 + k) + 2a/d] + rr2Wbhe2[t(2 + k/3) + 2a/d] 
+ re - 3d2b2R2[4(2 + 4/3) + 2a/d] = 0 


(10) 
Eliminating A; and Ag from the above 


bkt[3t(2 + k/3) + 4a/d] d2b2h2(2 + k/3) 2d2(2 + k) 
- bkt®[t(2 + k/3) + 3a/d] adb?h? 2ad 
2Qbkt2[t(2 + k/3) + 2afd] 3d2b2k2[t(2 + k/3) + 2a/d] 2d*(¢(2 + k) + 2a/d] 


Expanding 


IQh2 3 2 2 
aar(16 + 40H) (15 4 44 0), (5 28) a2 1008 0 


Putting td = ap the above equation may be expressed 
as :— 
aop*® + ap? — aep — az = 0 ~ (8%) 


Where 

@o = 72 + 60k + 14k? + FS 

@ = 72 + 102k + 15k? 

ag = 72 — 42k 

a3 => 72 
The two conditions (6) for the I-Section may be rewritten 
as :— 

(t/d)*b*k[t(2 + k/3) + 2a/d]) = M/2K’E - (12) 


b3k2[4(2 + k/3) + 2a/d]) = M/2E® ~~ ~- (13) 
Dividing (12) by (13) 


(¢/d)2 - PS = 0/K’ 


Substituting for ¢ and } in equation (13) gives 
b = 1-134[K’E — 4a— 20-52-13 g — 462M 41/18 
= 1-134[K’E — 4a — 26 — 5f,(k)M4)1/38 
= [K’'E — 4a — 26 —5M,4)1/18 
(14) 
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Where 
a = p + 6p/k + 6/k 
Ailk) = k-- 8a — 
M,=1- “13atiyerh ay eM 
Similarly solving for ¢ and d 
t = 1-032[K’ — 3E —1@% 62a — 197M}}/18 
= 1-032[K’ — 8E — 146 62/2(k)M)1/18 
= [K’ — 8E — 14% 62M9)!/18 


(15) 
Where 
fa(k) = a— 1p? 
Mz = 1-03213f2(k)M 
and 
d = 0-969[K’3 Ea? 0- 2ap6M —1)1/18 


Dividing 6 by d gives n. 


Hence 
n = 1-170[K’— 2E-5a— 9 6-3-1384 —5p- 4 5)1/18 
= 1-170[K’ — 2E —5a— 96 — 3f3(k)M5)1/18 
= (K’- 2E-5a— 99 — 34qs]I/18 
. . (17) 
Where 


fa(k) = k-18a—5p- 4 
Mz = 1-17018/5[ fg(k)}4/5M 
Channel Section 
By a similar reasoning the equation corresponding 
to 11 in the case of ‘channel section can be shown to be 
ao g° + aig? — aog— az = 0 . - (18) 
Where 
ao = 18 + 30k + 14hk2 + 2k8 
18 + 51k + 15k? 


a, = 
a2 = 18 — 21k 
a3 = 18 
= td/a 
and 
b = 1-134[K’E — 4a — 26 — 5k — 18g — 4924 4)1/18 
= 1-134[K’E - 4a- 29 —5g,(k) MS)}/13 
= [K’E- 4a- 26-5M,4}1/18 
(19) 
t = 1-032[K’ — 3E — 146 §26 —197M)1/18 
= 1 -032[K’ — 3E — 14% §2go(k)M)1/13 
= [K’— E -146 §2M9)1/18 
; (20) 
n = 1-170[K’-2E—5a— 9 §—3k—13g —5g— 4451/18 
= 1-170[K’ — 2E —5a- 96 — 3gg(k)M5}1/18 
= [K’— 2E-—5q— 99 — 3, 5)1/18 
(21) 
Where 
B=¢ + 3g/k + 3/k 
gi(k) = k—-136 - 
M,; = 1-134 19/4h5, (3) OM 
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g2(k) = B—1g? 

Mz = 1-032}3go(k)M 
&3(k) = k— 188 — 5g — 4 

Msg = 1-17013/5[g3(k)}}/5M 
Box Section 
The cubic equation in this case is 

aor® + ayr2 —agr—ag=0 - . (22) 

Where 


ao = 36 + 24k + 4h? 
@, = 36 — 15k — Qk? 


= 36 + 39k 
a3 = 36 
r = tdla 
and 
b = 0:915[K’E — 4a — 26 — 5k — 18y — 47274)1/18 
=0- pad E — 4a-2 6 — 5h, (k) M4)}1/18 
= [K’E — 4a-2 6 — 5M, 4}1/18 
t = 0-978[K’ — 8E — 146 62+ — 147gM}}1/38 
= 0-978[K’ — 3E — 148 62he(k)M}4/38 
= [K’— 8E —1q6 §2Mo}/18 
n = 0-895[K’—2E-5a— 9 @— 3h—-18y —5y— 4M5}1/18 
= 0-895[K’ — 2E —5a— 96 — 3hs(k)M5)1/18 
= [K’— 2E —5qa— 96 — 3M5}1/18 
Where 


=r-+ 3r/k + 3/k 
h s() = k—-18y — 4y2 
My = 0-91513/4(/4(&)}1/4M 
ho(k) = y— 1? 
M2 = 0-97813h2(k)M 
hg(k) = k- 187-5 — 4 
Ms = 0-89513/5[hg(k)]}/5M 


r 
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Discussion 


The Council would be glad to consider the publication of 
correspondence in connection with the above paper. Communi- 
cations on this subject intended for publication should be 
forwarded to reach the Instituticn by August 31st, 1961. 
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Theory and Practice of Arc Welding, 2nd Edition, by 
R. J. Sacks. (New York & London: D. van Nostrand 
Co., Inc., 1960). 9}in. x 6}in., 478 + xii pp. 45s. 

This is the second edition of a textbook first 
published in 1943, and is based on an original course of 
study in arc welding theory and practice in accordance 
with the current practices in industry in America and 
the recommendations of the American Welding Society. 

The first section presents fundamental welding theory 
and includes chapters on industrial use of arc welding, 


Book Review 






equipment, mild steel electrodes, joints and welds, the 
nature of iron and steel, testing and inspection, flame 
cutting and safety precautions. Questions are given 
at theend ofeachchapter. The second section presents 


a comprehensive welding training course and deals with 
manipulative skills. 

The text of the new edition has been revised to bring 
it in line with current practice, many new photographs 
and diagrams have been included and two new chapters 
have been added. 











May, 1961 


Institution Notices and Proceedings 


ORDINARY GENERAL MEETING 


An Ordinary General Meeting of the Institution of 
Structural Engineers was held at 11, Upper Belgrave 
Street, London, S.W.1, on Thursday, 23rd March, 
1961, at 5.55p.m. Lt.-Colonel G. W. Kirkland, 
M.B.E.(Mil.), M.I.Struct.E., M.1.C.E., in the Chair. 

The following members were elected in accordance 
with the Bye-Laws. Will members kindly note that the 
elections, as tabulated below, should be referred to 
when consulting the Year Book for evidence of 
membership. 

STUDENTS 
BROWNE, Geoffrey, of Luton, Beds. 
CocKBILL, David, of Salisbury, Southern Rhodesia. 
FERNANDO, Kalutarage Stanley, of Colombo, Ceylon. 
Hype, Stephen Arthur, of London. 
Jounson, Alan Edward, of Wilmslow, Cheshire. 
MBAEZUE, Joseph Ikecukwu, of Zaria, Nigeria. 
Parks, Anthony, of Middlesbrough, Yorks. 
SatIM, Mohammad, of District Sheikhupura, West 

Pakistan. 

SRIVASTAVA, Upendra Nath, B.Sc., of Gorakhpur, 
- Uttar Pradesh, India. 

WaLKER, David Whinray, of Morecambe, Lancs. 

Weo, Peter, of Middlesbrough, Yorks. 

Woop, Jonathan Guy Martin, of Frensham, Surrey. 


GRADUATES 

ApaMs, James Edward, of Sutton, Surrey. 

ARTHERN, Frederick, of Stockport, Cheshire. 

BEVAN, David Milton, of Winterbourne, Bristol. 

Brown, Brian Arthur, of Winchester, Hants. 

CHANDRACHOOD, Madhukar Abaji, B.E., of Poona, 
India. 

CoLES, Brian Stanley, B.Sc., of Horley, Surrey. 

CROCKER, John, of Loughborough, Leics. 

Desal, Chandrakant Sankalchand, B.E., of Bombay, 
India. 

Duar, Chaman Lal, B.Sc., of Bombay, India. 

FroGGATT, Michael Trevor, of Leamington Spa, 
Warwicks. 

GranBy, Alan Frederick, of Liverpool, Lancs. 

GREEN, Michael Henry, of Kingsworthy, Nr. Win- 
chester, Hants. 

HENNEY, Graham, of West Wickham, Kent. 

Hor_er, David Brian, of Whitchurch, Cardiff. 

ImEGwu, Eugene Odili, Ph.D., B.Sc., of Bristol. 

Japav, Manhar Kalidas, B.E., of London. 

Kone Func Kew, M.Sc., of Leeds. 

MarTIN, Thomas Allan, B.Sc., of London. 

Mitra, Gopal Chandra, B.Tech.(Hons.), of Cologne, 
Western Germany. 

Nig, Peter, of London. 

O’CALLAGHAN, Michael John, B.E., of Newport, Mon. 

Rance, Ronald John, of Tring, Herts. 

SHER Kar Kit, Joseph, of Holmcroft, Staffs. 

SOLANKI, Chhaganlal Dewji, B.E., of London. 

STEER, Peter John, B.Sc., of Hounslow, Middlesex. 

Tuan, Aung, B.Sc., of Rangoon, Burma. 

THIRUMAL, Paramalingam, B.Sc., of Kopay, Ceylon. 

Totat, Vipin Natwarlal, B.E., of London. 

Upgaja, Agbasielo Pius, of London. 


ASSOCIATE-MEMBER 
Rowe, Roy Ernest, M.A., A.M.I.C.E., of Slough, Bucks. 


MEMBER 
McNutty, Jim, of East Didsbury, Manchester, 20. 


TRANSFERS 
Students to Graduates 

Basu, Utpal Kumar, B.Sc., of Greenford, Middlesex. 
BisHop, John Leighton, of Solihull, Warwickshire. 
CLEATOR, Michael John, of Harrow, Middlesex. 
FULLER, Norman, of Manchester. 
Goopwin, David Laurence, of Sutton, Surrey. 
HAWKER, Brian Howard, of Cardiff, Glam. 
Hooper, Jeremy, of Ilkley, Yorks. 


Graduates to Associate-Members 
ASKEY, Reginald Charles, of London. 


a 


Associate-Members to Members 

BrIMER, Alexander, of Lagos, Nigeria. 

BuLtock, Dennis Terence, B.Sc., A.M.I.C.E., of Sutton 
Coldfield, Warwicks. 

BuncE, John Wallace, of East Ewell, Surrey. 

CuRSITER, William Nisbet, B.Sc., A.M.I.C.E., of 
Lagos, Nigeria. 

Davis, Frank, of Mombasa, Kenya. 

FRAENKEL, Peter Maurice, B.Sc.(Eng.), of London. 

HARVEY, John Alan, B.Sc., A.M.I.C.E., of Salisbury, 
Wilts. 

HERRING, Bernard Carroll, of Lagos, Nigeria. 

McNAUGHTON, Finlay Malcolm, A.M.I.C.E., 
A.M.I.Mun.E., of Bolton, Lancs. 

Moran, Thomas Francis, of Liverpool. 

PINKNEY, John Dennis, of Scunthorpe, Lincs. 

VERYARD, Joseph Gordon, of Gabalfa, Cardiff. 


RE-ADMISSION 
Assoctate-M embership 
GREGORIAN, Goriune Thaddeus, of Sidcup, Kent. 


OBITUARY 
The Council regret to announce the deaths of Charles 
Dimond Conway BRAINE (Member), and Leslie Guest 
Harris (Retired Associate-Member). 


RESIGNATIONS 

Notification was given that the Council had accepted 
with regret the resignations of Raymond Victor CHATE 
(Member), William McKay (Jnr.) (Associate-Member), 
Saghayroon Et ZEIN, William Thexton LaAwTOoN, 
Bent Mehl LupvicsEN, Ebenezer Olufunso ODEDAIRO, 
Michael John Puncu (Graduates), Thomas Neville 
WELLs (Student). 


FORTHCOMING MEETINGS 


The following meetings will be held at 11, Upper 
Belgrave Street, London, S.W.1. 


Thursday, 25th May, 1961 
Ordinary General Meeting for the election of 
members at 5.55 p.m., followed by the Annual General 
Meeting at 6 p.m. 
Thursday, 22nd June, 1961 
Ordinary General Meeting for the election of 
members 5 p.m. 


INSTITUTION GF STRUCTURAL ENGINEERS 
BENEVOLENT FUND 
The Annual General Meeting of the Voting Contri- 
butors to the Institution of Structural Engineers 
Benevolent Fund will be held at 11, Upper Belgrave 
Street, London, S.W.1, on Thursday, 25th May, 1961, 
at 6.20 p.m. 





EXAMINATIONS, JANUARY, 1961 


HOME AND OVERSEAS CENTRES 


The Examinations of the Institution were held in 
January, 1961, at the usual centres in the United, 
Kingdom, and overseas at the following centres :— _ 

Accra, Aden, Aligarh, Auckland, Beirut, Bengazi, 
Bombay, Brazil, Cairns, Cairo, Calcutta, Cape Town, 
(South Africa), Christchurch (New Zealand), Colombo, 
Dar-Es-Salaam, Durban (South Africa), Harrisburg, 
Hong Kong, Johannesburg (South Africa), Kampala, 
Karachi, Khartoum, Kingston, Komati, Kuala Lumpur, 
Lagos, Lucknow, Lusaka, Madras, Mauritius, 
Melbourne, Montreal, Nairobi, Peru, Port Elizabeth 
(South Africa), Salisbury (Southern Rhodesia), 
Singapore, Sydney, Toronto, Trinidad, Vancouver, 
Wellington, (New Zealand) and Windhoek (South Africa). 

One hundred and sixty-two candidates took the 
Graduateship Examination (103 at home and 59 
overseas), and 808 candidates took the Associate- 
Membership Examination (650 at home and 158 
overseas). Of these, 81 passed the Graduateship 
Examination (44 at home and 37 overseas) and 154 
passed the Associate-Membership Examination (135 
at home and 19 overseas). 

The names of the successful candidates are : 


Graduateship Examination 


ApEsoca, Mobolaji 

AFOLABI, Joseph Oduntan 
ARTHUR, Roland Eric Brough 
ATHALYE, Shridhar Mohaniraj 
Buat, Kaje Rama 

CaTo, Ross Edward 
CHABLANI, Mohan Ghanshamdas 
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Cu1u CHI SHING 
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FARNWORTH, Frank 
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FERNANDEZ, Osborne Anthony 
ForrEsT, Frank D. 
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GLOVER, Benjamin Payne 
GODFREY, Dan 
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GRIFFITHS, John Charles 
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HAQuE, Malik Israrul 
HinGorANI, Kishin Bhagchand 
HuGueEs, Peter James 
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ISANEDIGHI, Okon Uwe 

Lal SHUI KWAN 

LARNER, Raymond James 
LEE CHEOK HOE 

LockE, Joseph 

LOREILHE, Vasco 

McCarTAN, Peter John 
MaGuIRr_E, Robert John 
MANIAR, Priyavadan Umedlal 
MILLER, Roy James 

Moores, Arthur Keith 
NARASINGA Rao, Prodduturi 
Notice, Ezekiel Stephen 
Os1ozo, Leo Emmanuel 
OcGuNSOLA, Thompson Olusola 
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OKAFoR, Ebenezer Chikwudozie 
OLIVEIRA, Denis Ernest Ligorio 
OmER, Awadalla Mohamed 
Onvouna, Nicholas Onyemauna 
OweEN, John Alan 

Pao WEI KANG 

Parsons, Eric Henry 

PATEL, Kantilal Purushottamdar 
PAUL, Iqbal Singh 

Pieris, Kurukulasuriya Justin Joseph S. 
Pitot, Miss M. A. G. 

POTGIETER, Theodorus Daniel 
RaHA, Jogalrata 

REID, Michael Irwin 

Rizvi, Aulad Hussain 
ROTHWELL, Conrad 

SAFRI, T. L. 

SALAKO, James Adebayo 
SCHROEDER, Werner Otto Esmarch 
SCHULER, Ronald Manfred 
SESHAGIRI Rao, Kompella Venkatu S. 
SHARMA, Har Bhagwan 

SHASTRI, M. K. 

S1ncH, Dharam 

SIVALINGAM, Gunaretnam 
SmiTH, Peter William 

TAanG Hunc CHUN 

Tuomas, Warwick 

TOwWLER, Malcolm Philipson 
VAMADEVAN, Nagalingam 

WEBB, Ernest George 

WESSIER, Bernard Eugene 
WHELAN, Peter Denis 

Wone CHING WAH 

Wricut, Brian Allen 

WRroE, Herbert Brian 

YunG, Oliver Hon Bun 


Associate-Membership Examination 


ABu-SITTA, Salman Husein 
ANDERSON, David Graham 
Avis, Peter Derek 
BEAMAN, Jack 
BENNETT, John David 
BHOoGAL, Harbans Singh 
BisHop, Colin Charles 
BLASszyYK, Tadeusz 
Boyce, Michael David 
BricuT, Norman John 
Brown, Roy Charles 
BULL, Geoffrey Michael 
Bunn, David Michael 
BurRNETT, Walter 

Busu, Kenneth Donald 
CARSON, Paul Frederick 
CHANDA, Kanak Kumar 
CHATARJEA, Sukhenda 
CuHoprRA, Narindar Nath 
CLARKE, Francis Gerard 
CLEMENTS, Alfred Cecil 
CLoucu, Norman 
CORBETT, Brian Oliver 
Corry, John Percival 
Cox, Anthony James 
CoyYNeE, Terence 

Cross, John Henry 
Cruppas, Edward Alan 
DEAKIN, William Herbert 
Dixon, Raymond James 
Dopp, Harry Thomas 
Dopp, James Michael 
DovEtT, Louis Maurice 
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Dovc.as, James P. D. 
DrakKE, Michael John 

Dutt, Nihar Kusum 

EaALEY, Michael John 
Epwarps, William John Stephen 
Ex.ison, Morville Vernon 
Evans, David Geoffrey Rhys 
FAIRWEATHER, Roy Leonard 
FARSIDES, Costa Pavlou 
FitzjoHN, Allan Arthur 
FLANAGAN, George Edgar 
FLEMING, Ian Doyle 

FRASER, Frank Kojo 
FRIBBANCE, William Lawrence 
GacH, Jerzy 

GHosH, Srikumar 
GILLIBRAND, William 

Gray, Bruce Alan 

GREAVES, Brian Eric 
GRIEVE, Neil Macpherson 
GROVE, James Arthur 
Happock, Geoffrey 
HAGENBACH, Paul 
HALESTRAP, Terence Leonard 
HALL, Conrad 

HALLET, Rex 

Harris, Clifford James 
Harris, Martin Anthony 
Harris, Peter William 
HasTINGS; Joseph William 
HERON, James Alexander Stewart 
HojEM, Peter Anthon 
HoLroyD, Philip 

Hottum, Ronald David 
HookINs, Barrie Scott 
HuFFORD, Howard Horace 
HuGHEs, John Allan 

Hunt, David Charles Gowar 
HunTeER, Ian Henderson 
James, Harold Cecil 

Jay, Bernard Ensor 

Jerson, Wilfrid John 
JorDAN, Edward John 
Kerns, Donald Patrick 
KinG, Peter Henry 
Kulkarni, H. M. 

Lainc, John Cranston 
LAWRENCE, Hedley 

Lawson, Harland 

Laxton, David James 

Li Kur Wal 

LEIBNER, Rolf Otto 

LIMNER, Thomas 

Lipp, Lawrence 

LockE, Joseph 

LoucuTton, Alan James 
LovELEss, Colin 

LupLaM, Donald Thomas 
Lunar, J. D. 

McGratTH, Frederick James 
McKay, Colin 

McLean, Robert 

MacRag, Herbert Tweddle 
Manin, Carl Lionel 

Manayan, K. D. 

MARKLAND, Frederick George 
MaTHEw, George 

Micuna, Tadeusz 
MIDDLEMISS, Joseph Watson 
MILLIGAN, Owen Page 
Munscu, Anthony Mark 
Nice ILLS, Tose np} Thomas 


NIGHTINGALE, Peter 
Nurse, Colin Frederick 
NurseEy, John Richard 
ORENSTAIN, Jehuda Leo 
OXLEY, Victor John 
PALMER, David Ernest 
PATTERSON, Robert Alfred 
PERINI, Flavio 

PiLGRIM, Ernest Joseph 
PITCHERS, John 

PonTIN, Ronald Frederick 
POOLE, James John 
PRABHU, Sudhakar Shrivang 
Rawson, Morris 
RICHARDSON, Albert 
RICKETTs, Daniel John 
RosBInson, Frank Leonard 
Rosinson, George Leslie 
ROBINSON, Michael Douglas 
SANDILANDS, Brian Sidney 
SATHE, R. M. 

SEN GupTA, Pranab Kumar 
SETT, Deb Kumar 

SLADE, Harry Jerrold 
SMITH, John Frank McLaren 
SPICER, Douglas Aubrey 
STIEFEL, John Emil 

STONE, Peter Albert 
SyMES, Ronald Clive 
TAYLOR, Eric Henry Frank 
THomsoN, John Stuart 
TosutTt, David Cooper 
Trump, Vernon Henry 
TURNBULL, James Neil 
TURNER, Frank Herbert 
WALKER, Derek 

WALLACE, Reginald Rutherford 
WEISE, Milton Ludlow 
WELsH, John Sturgeon 
WEsTON, Paul Anchorage 
Wuite, Keith Christopher 
WuiTE, Stanley Francis 
WuitinG, David Stanley 
WILImMowskKI, Zbigniew Marian 
WItson, Frederick 
Witson, Norman Angus 
Woop, Robert Francis 
Younc, James 

ZEEGEN, Alan Stanley 


EXAMINATIONS, JULY, 1961 


The Examinations of the Institution will next be 
held in the United Kingdom and Overseas on the 11th 
and 12th July, 1961 (Graduateship), and the 13th and 
14th July, 1961 (Associate-Membership). 


CHANGE OF TELEPHONE NUMBER 


The telephone number of the Institution has been 
changed to BELgravia 4535 (3 lines). 


YEAR BOOK AND LIST OF MEMBERS 


The Year Book and List of Members for 1961 will 
go to press in July, for publication in October, when 
a copy will be sent to all members. 

Members are requested to inform the Secretary of 
any alterations in titles, degrees or addresses, which 
have not already been notified, by June 24th, in order 
that such amendments may be included in the new 


edition 
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COMMONWEALTH TECHNICAL TRAINING 
WEEK 


The Institution is staging an exhibition of photo- 
graphs, models and drawings illustrating the profession 
of Structural Engineering from 29th May to 2nd June, 
at 11 Upper Belgrave Street, London, S.W.1. 

The exhibition will be open daily from 10 a.m. to 
12 noon and from 2p.m. to 8p.m. Members of the 
Institution will be available throughout to answer 
questions regarding training for the profession and 
prospects. 

During the mormings of Monday, 29th May, 
Wednesday, 31st May, and Friday, 2nd June, coaches 
will be provided to take those interested to visit 
outstanding structures under construction such as 
Hammersmith Flyover, the Royal College of Arts, 
and some of the new high buildings. 

During the afternoons of these days illustrated 
talks will be given on various aspects of Structural 
Engineering and on the Tuesday and Thursday there 
will be film shows. 

It is hoped that members and their friends who are 
able to do so will visit the exhibition, particularly 
between 6 and 8 o’clock each evening. 


DRURY MEDAL AWARD 


The alternative subjects for this eighth competition 
are a pipe-bridge’and a petrol service station. 

Graduates and Students of the Institution not over 
27 years of age are invited to apply for full details 
to the Secretary, envelope to be marked at the top 
left-hand corner ‘ Drury Medal Award.’ 

The general conditions of the competition are as 
follows :— 


(a) The competition shall be for Graduates and 
Students of not more than 27 years of age. 


(b) The subjects of the competition will be designs 
of a structural character, that is to say, involving 
structural design rather than planning. 


The subjects of design and the conditions shall 
be prepared and issued biennially. 


A Jury shall be appointed to examine the work 
submitted and to interview candidates if found 
necessary. 


In order to ensure that the design submitted 
is the unaided work of the competitor, the 
drawings, calculations, etc., submitted shall be 
endorsed by the candidate: “I declare that 
the work I hereby submit is my own unaided 
work.” The declaration shall be signed by the 
competitor, and be either countersigned by a 
corporate member, or be certified as made 
before a Justice of the Peace, or a Commissioner 
for Oaths. 


Branch Notices 


LANCASHIRE AND CHESHIRE BRANCH 
MERSEYSIDE PANEL 


The Annual Dinner-Dance will be held at the Carlton 

Chub, Eberle Street, Liverpool, on Friday, 12th May, 
1961. 
Joint Hon. Secretaries: Wm. Watts, M.LStruct.E., 
A.M.LC.E., 11, Newchurch Lane, Culcheth, Nr. 
Warrington, Lancs., and M. D. Woods, A.M.I.Struct.E., 
8, Dennison Road, Cheadle Hulme, Cheshire. 


The Structural Engi . 


MIDLAND COUNTIES BRANCH 


Hon. Secretary : S. M. Cooper, M.I.Struct.E., ‘ Ap 

garth,’ 56, Hyperion Road, Stourton, Stourbridgs 

Worcs. Bi 
GRADUATES’ AND STUDENTS’ SECTION . 


Hon. Secretary: H. T. Dodd, Shepherd’s Cottag 
Grove Lane, Wishaw, Sutton Coldfield, Warwicks, © 


NORTHERN COUNTIES BRANCH 


Hon. Secretary: P. D. Newton, B.Se.(En ), 
M.I.Struct.E., A.M.I.C.E., 6, Cornfield Road, Linthorp 
Middlesbrough, Yorks. 


NORTHERN IRELAND BRANCH 


Hon. Secretary: L. Clements, A.M.1.Struct Ey 
A.M.I.C.E., A.M.I.Mun.E., 3, Kingswood Park, Cher 
valley, Belfast, 5. i 


SCOTTISH BRANCH 


Hon. Secretary: W. Shearer Smith, M.1StructE 
A.M.I.C.E., c/o The Royal College of Science ap 
Technology, George Street, Glasgow, C.1. 


SOUTHERN BRANCH a 

Hon. Secretary: A. P. K. Tate, B.Sc., A.M.I.Struct.E 

Department of Civil Engineering, The University 

Southampton. 4 
SOUTH-WESTERN COUNTIES SECTION 


Hon. Secretary: C. J. Woodrow, J.P., ‘ Elste . 
Hartley Park Villas, Mannamead, Plymouth, De 


WALES AND MONMOUTHSHIRE BRANCH 


Hon. Secretary: W. D. Hollyman, A.M.I.Stru 
41, Greenfield Avenue, Dinas Powis, Glam. 


WESTERN COUNTIES BRANCH 


Hon. Secretary: A. C. Hughes, M.Eng., M.L.Stru 
A.M.I.C.E., 21, Great Brockeridge, Bristol, 9. 


YORKSHIRE BRANCH 


Hon. Secretary: W. B. Stock, A.M.I.Struct.E., n 
Hobart Road, Dewsbury, Yorks. a 


UNION OF SOUTH AFRICA BRANCH ~ 
Hon. Secretary: E. B. Kretzchmar, A.M.1.Structl 
P.O. Box 3306, Johannesburg, South Africa. 
Natal Section Hon. Secretary: J. C. Pant 
A.M.1.Struct.E., A.M.I.C.E., c/o Dorman Long (Affe 
Ltd., P.O. Box 932, Durban. 7 
Cape Section Hon. Secretary: R. F. Nom 
A.M.1.Struct.E., African Guarantee Building, 8, 3 
George’s Street, Cape Town. 4 


EAST AFRICA SECTION 


Hon. Secretary: K. C. Davey, A.M.I.Struct.E., 
Box 30079, Nairobi, Kenya. 


SINGAPORE AND FEDERATION 
OF MALAYA SECTION ty 
Hon. Secretary: J. R. M. MacIntyre, A.M.1.Struet 
c/o Redpath, Brown & Co. Ltd., P.O. Box 
Singapore. 
NIGERIAN SECTION 
Hon. Secretary: A. Brimer, M.I.Struct.E., Br 
Andrews & Nachshen, Private Mail Bag 2295, 
Nigeria. 3 
AUCKLAND (NEW ZEALAND) SECTION © 
Hon. Secretary: A. Donald, 3B.Sc.(Hom 
A.M.I.Struct.E., 122, Matipo Road, Te 
Auckland, New Zealand. 





ae 





